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SYNOPSIS 

The present investigation deals with the study of the 
second-order laminar boundary- layer equations. The three 
aspects of the problem studied are; (i) general self-similar 
solutions including dissipation,, (ii) locally-similar solutions 
near a stagnation point using Gortler method with an application 
to a blunted wedge and (iii) locally-similar solutions valid at 
large distances in the downstream direction with an application 
to a blimted wedge. 

In the first problem, we investigate a general structure 
of the self-similarity for the second-order boundary-layer 
equations with dissipation. The second-order boundary-layer 
equations associated with longitudinal curvature, transverse 
curvature, displacement speed, external vorticity and enthalpy 
gradient are obtained from the Navier-Stokes equations by using 
the method of matched asymptotic expansions. As regards the 
thermal boundary condition on the wall we study both cases; 



xx 


that of the prescribed wall temperature and of the insulated 
wall. For the latter case, the classical concept of the 
recovery factor is extended to the second-order flows. The 
self-similar governing equations are integrated numerically 
and results are critically discussed; in particular, the 
present work explains the mathematical as well as physical 
aspects of the singularities observed in the second-order 
boundary-layer equations (See Afzal and Oberai 1972, Werle 
and -Davis 1970). In the foimulation of the thermal problem, 
the thei?mal expansion parameter |~defined as (aiogf/aiogT)^, 
where f is the density, T is the temperature and p is the 
prsssur^ is negligibly small and hence the comp.ressibility 
term is neglected. This is valid for most of the liquids. 
However, for the flow of perfect gases, the magnitude of the 
thermal expansion parameter is unity and thus the compressibility 
term is to be retained. For such cases, a general formulation 
of the second-order boundary-layer equations is given by 
Van Dyke (1962), Self-similar solutions to Van Dyke's equa- 
tions are also obtained. 

Once the general structure of the self-sinilar solutions 
is understood, the next logical step is to develop nonsimilar 
solutions. The second problem deals with locally- similar 
solutions for a unifomn flow past two-dimensional bodies -wtoere 
only displacement and longitudinal curvature effects are present. 
The method is illustrated by solving the problem of a blunted 
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wedge IjLf denotes the wedge angle then parabola and flat 

plate normal to flow direction are special cases of the blunted 
wedge with = 0 and 1' respectively]. The displacement problem 
is global in the sense that it alone req.uires the determination 
of the outer flow due to displacement. The perturbed Euler 
equation is solved and the displacement speed, U 2 (<t '5 O) , is 
obtained as a function of streamwise coordinate , <cr' . First five 
terms for the se'cond-order skin friction and heat transfer in a 
Gdrtler series expansion in cv are calculated. The convergence 
of the series is poor for large cy • The region of validity is 
extended in the downstream direction by Eulerization of the 
series. Non-linear Shanks transformation has also been employed, 
where necessary, for further accelerating the convergence. It is 
found that the Eulerized results for the skin friction and heat 
transfer give accurate results (within 2 percent) even as cr* °o. 
Further, a comparison of the present results for the special 
case of a parabola, with the exact solutions of the Favier- 
Stokes equations of Davis (1972), Dennis and Walsh (1971) and 
Botta, Dijkstra and Veldman (1972) show a good agreement (within 
one percent) upto Reynolds number as low as one thousand. 

The series solutions of the kind mentioned above have, in 
general, a restricted region of validity in downstream direction 
and, therefore, it seems desirable to have alternate solutions 
valid for large y ; this forms the subject matter of the third 
problem. Solutions are constructed for the second-order 
boundary- layer equations for cr ~ by expanding the variables 
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in terms of its eigenfunctions with certain unspecified multi- 
plicative coefficients, to be determined from the prescribed 
initial conditions. The eigenvalues and eigenfunctions for 
various second-order effects are obtained by solving the asso- 
ciated homogenous problems, the governing equations of which 
have a foim similar to the well-known Sturm-liouville equation. 

Por the first-order bo undaiy- layer equations, such studies 
have been made by Libby and his coworkers (refer Libby 1970). 

The solutions for the second-order effects are displayed graphi- 
cally and discussed for a variety of flow situations. The results 
provide a further insight into the singularities of the second- 
order solutions. The results of the first problem explain the 
location of isolated singularities whereas those of the last 
problem show two distinct domains of stable and vinstable solu- 
tions separated by the first critical value (corresponding to 
the above mentioned singularities for various second-order 
effects) of the parameter. Thus, it is shown that if the second- 
order parameters assume a value less than the critical value, 
any departure from strict similarity grows in the downstream 
direction and as such the flow is, in some sense, unstable. The 
application of the above procedure is illustrated by considering 
a uniform flow past a blunted wedge. The initial condition is 
prescribed in terms of the Eulerized series solution obtained 
earlier (second problem). In addition to eigenfunctions, we 
include terms due to variation of the second-order parameters 
with C 3 -. The multiplicative coefficient of the first eigenfunction 
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IS determined for each of the second-order effects, by joining 
the present series - .00) with the Eulerized Gortler series. 
The resulting asymptotic solution for large .t' , in conjunction 
with the Eulerized Gortler series may be used to assess the 
accuracy of the corresponding results obtained by numerical or 
approximate methods, where the error bounds are unkno?/n or 
uncertain. 



CHAPTER 1 


INTRODCCTION 

Recently a lot of interest has -been shown in the 
so-called higher order approximations of the houndary-layer 
theory (specially the second-order). There are various 
reasons for undertaking such studies. Por example, the lead- 
ing term in the asymptotic expansion of the Havier-Stokes 
eq,uations for large Reynolds number represents the classical 
boundary-layer equations of Prandtl. It is well-known that 
the Navier-Stokes equations are elliptic whereas the classical 
boundary-layer equations are parabolic and thereby do not 
permit the upstream influence. It is of interest to study 
how the higher order terms in an asymptotic expansion reassert 
the elliptic nature of the Navier-Stokes equations, suppressed 
in the classical boundary-layer equations. The practical 
utility of such higher order corrections lies in the fact that 
they are found to give good results even under distinctly 
nonasymptotic situations; Lagerstrom and Cole (1955) have shown 
that the second approximation in some cases may predict correct 
skin friction upto Reynolds number as low as ten. 

There is another motivation for studying the higher 
order approximations. Prandtl 's boiandary- layer theory does 
not account for what Rott and lenaird (1959) call the secondary 
effects. Physically these effects can arise due to longitudinal 
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ciirvature, transverse curvature, displacement speed, external 
vorticity, external enthalpy gradient, velocity slip and tem- 
perature jump. These effects are of order as compared 

to the Prandtl's houndary-layer and become significant at 
moderately large values of Reynolds number. 

Van l^rlce ( 1962) gave a systematic formulation of the 
scheme for calculating the higher order corrections to the 
classical boundary-layer (first-order)* The technique used is 
that of matched asymptotic expansions developed by Kaplun, 
Lagerstrom and their colleagues* This method divides the 
flow into inner and outer regions which are to be matched in 
an overlapping domain. An interesting feature of perturbation 
technique is that the resulting second-order equations are 
linear and, thus, the analysis can be further simplified by 
subdividing them into various simpler problems, each having 
a clear physical interpretation. 

The theory so formulated by Van I^rke ( 1962) has been 
studied by various authors and an excellent review has been 
given by Van lyke (1969) and as such will not be repeated here. 
Among the articles published thereafter, the following three 
are most relevant for the present study: Afzal (1969), Werle 
and Davis (1970) and Afzal and Oberai (1972). These authors 
have proposed a more general form of the self-similar solutions 
for the second-order equations of Van lyke by identifying three 
new parameters , one each for the longitudinal ctjrvature , transverse 
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ciqrvature.and displacement speed. Werle and Ilavis(l969) studied the 
problem of momentum transfer while Afzal (1969) and Afzal and 
Oberai (1972), have also included the study of the heat transfer 
but without the viscous dissipation terms. These authors re- 
ported that for some critical values of the second-order 
parameters, the second-order solutions show several singulari- 
ties (infinite discontinuities), but the general pattern of 
the singularities was not investigated. The present work 
(Chapter 3) elaborates the mathematical as well as physical 
aspects of such singularities. It is shown that the singulari- 
ties occur whenever the parameters take the value equal to the 
eigenvalue for the corresponding homogeneous problem. 

The problem of heat transfer from two-dimensional and 
axisymmetric geometries has considerable practical interest. 

Afzal and Oberai (1972) in their study of heat transfer have 
neglected the viscous dissipation. However, there are situa- 
tions where the dissipation terms may not be negligible e.g. 
the case where the product of the Prandtl nixmber <T and Eckert 
number E is of the order unity i.e. or* B '^0(1). One such 
situation arises with liquids of moderately large Prandtl 
number (of the order of 50 or so). The formulation of Van 
Dyke (1962) for the thermal problem is applicable to perfect 
gases as the formulation takes care of the equation of state , 
i.e. P =:pET where P is the pressure, P the density, 1 the gas 
constant and T the temperature. Thus, for the flow of liquids, 
a new formulation is needed and is described in Chapter 2. 
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One of the problems studied in the present work deals 
with the second-order effects on heat transfer including 
dissipation for both liq.uids (employing the formulation of 
Chapter 2) and gases with constant density (using the formula- 
tion of Van Dyke). As regards the the3:mal boundary condition 
on the wall we study both cases: that of the prescribed (arbitrary) 
wall temperature and of the adiabatic wall (no heat transfer at 
the wall). For the latter case, the classical concept of the 
recovery factor is extended to the second-order flows, A general 
self-similar analysis of the heat transfer with dissipation for 
both liquids and gases is presented in Chapter 3, The equations 
so obtained for various second-order effects are integrated 
numerically for various values of parameters by Runge-Kutta- 
Gill method with a step sizeA'^- 0.5 on IBM- 7044 Computer at 
IIT Kanpur. The results are displayed graphically and a criti- 
cal estimate of contributions due to various second-order effects, 
attributed to longitudinal curvature , transverse curvature , 
displacement speed, external vorticity and enthalpy gradients is 
presented. 

Once the general structure of self-similarity for the 
second-order boundary-layer equations has been understood, the 
next logical step is to develop the methods for studying the 
nonsimilar solutions. To deal with nonsimilar flows, various 
methods are available in literature, e,g, numerical techniques 
(like finite difference, finite elements etc,), momentum 
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integral methods and series solutions using some sort of local 
similarity. The numerical methods like finite difference 
employed for studying the full partial differential equations 
are expensive and need special experience. Further, in a given 
problem a fresh calculation is needed for each set of data so 
that the effects of the parameters on the solution are not 
obvious. The momentum integral methods are approximate methods 
and, in general, it is difficult to estimate the accuracy. The 
methods for studying the locally-similar flows (starting from 
the stagnation point or leading edge) have been constructed by 
Blassius (1908), Gortler (1957) and Meksyn (196l), The technique 
of Meksyn was further developed by Merk (1959) and Bush (1964) 
and will be cited as MMB method. It is well-known that the 
Blassius method suffers from the disadvantage of poor conver- 
gence of the series for a slender body shape (cf.e.g. Schlichting 
4968). The G-drtler and the MMB methods are attempts to overcome 
the above mentioned drawback of the Blassius method. In the 
Gbrtler and MB methods, the leading term (governed by the well- 
knovm Falkner-Skan equation) satisfies the outer boundary condi- 
tions exactly while the other terms of the locally-similar 
expansion give corrections that arise from the inner part of 
the boundary-layer. The difference in the G-6rtler and the MMB 
methods lies in the choice of the expansion parameter. In the 
Gdrtler method, the expansion parameter is 

s 

ex's: f U,(s,o) ds 
0 
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while in the MMB method, it is dp/djjjf wtoere 

p as 2d In U^/d In CT and X = Vo^ 

In the MMB method, the numerical solution involves terms like 
d^f^dpdX for r = 0,1,2. . .(here is the similar 
solution). Ihus, the calculations tend to he tedius and 
laborious. Considering all the aforementioned methods, it 
seems that a compromise between the accuracy and the computa- 
tional labour involved is achieved by GSrtler type series 
solution and is used here for the study of the nonsimilar 
flows considered next. 

In the present work we apply the above method to study 
the locally similar solutions of uniform flow past two-dimen- 
sional bodies. In such flows, the second-order effects are 
those of longitudinal curvature and displacement. Both momentum 
and heat transfer are studied. Using the G-drtler method, 
solutions are obtained by expanding the flow variables in terms 
of the streamwise coordinate, o'. It seems desirable that the 
solutions for the general case be obtained by constructing the 
universal functions which when multiplied by the appropriate 
parameters of a given problem, yield the required results. 
However, it was observed that such a procedure leads to unmana- 
geable number of equations and, therefore, no attempt has been 
made to tabulate the universal functions. Instead, the procedure 
is fully illustrated for a particular problem, namely flow past 
a blunted wedge. 
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We assume a steady laminar flow, at moderately large 
Reynolds number, past a semi- infinite, symmetric blunted wedge 
(shown in Figure 4.1) with constant wall temperature. If (Poo'^) 
be the included angle of the blunted wedge, two well~known 
special cases of the blunted wedge are parabola ( Poo = O) and 
flat plate noimal to flow direction ( - l). For a blimted 
wedge, the first-order momentum transfer has been studied by 
Chen, Iiibby, Rott and Van Dyke (1969) while the corresponding 
heat transfer is studied by Chen (1970). The second-order effects 
for a parabola (p^o = ®) been studied by Van Dyke (1964a) by 

estimating the leading term in the Blassius series expansion and 
as such the second-order contributions are evaluated at the 
stagnation point only. 

Of the two second-order effects, that due to the longitu- 
dinal curvature is a local effect as it depends only on the ■ 
first-order boundary- layer quantities while the displacement 
effect is global in the sense that it requires the evaluation 
of the outer velocity, U2 (s,o) and thus, represents the elliptic 
nature of the original Navier-Stokes equations. In the present 
work, we obtain the required displacement speed, T72(®>®)5 
solving the perturbed Euler equation with appropriate boundary 
conditions at far upstream and on the body surface. However, 
there is an alternate procedure to find the correction in the 
outer flow due to displacement thickness. The Euler equation 
is solved for a fictitious body with appropriate boimdary 
conditions. The profile of the fictitious body is obtained by 
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a suitable fit to the original body plus its displacement thick- 
ness (cf. e.g. Davis and Plugge-Dotz 1964). 

Once the expressions for displacement speed and longitu- 
dinal curvature are known, we obtain the governing equations 
for the first five terms of a series in cr* , for both skin fric- 
tion and heat transfer. These equations are integrated numeri- 
cally for various values of wedge angles. However, the conver- 
gence of the series is poor for large cr , This happens, often, 
due to the presence of a singularity in the complex-plane and 
as such the radius of convergence is limited due to a mathema- 
tical singularity rather than a physical one. In such cases, 
it is possible to enlarge the radius of convergence by either 
removing the singularity or pushing it farther downstream. One 
of the most powerful techniques for such purposes is Euler iza- 
tion which is employed when the location of the singularities 
is either known or can be guessed. Eurther, if there is only 
one singularity and its location is known, the Eulerization 
extends the radius of convergence to infinity. Vfe have foxuad 
the location of the singularity for our case so that the power 
series of skin friction and heat transfer is Eulerized to extend 
the validity range farther downstream. It is sometimes possible 
to further accelerate the convergence of a series by either 
applying the Shanks non-linear transformation or using extra- 
polation, These have been employed, where necessary. It is 
fotmd that the Eulerization, combined with extrapolation, gives 
accurate results (within two percent) for the skin friction and 
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heat transfer even as 0^«>. further, a comparison of the present 
results for the specis.1 case of a parabola (^^0 = 0), with the 
exact solutions of the Navier-Stokes equations of 3)avis (1972), 
Dennis and Walsh (1971) and Botta, Dijkstra and Veldman (1972). 
j^ows a good agreement (within one percent) upto Reynolds number 
as low as one thousand. 

The series solutions of the kind mentioned above have , 
in general, a restricted region of validity in the downstream 
direction and, therefore, it seems desirable to have alternate 
solutions valid for large <sr , This forms the subject matter 
of Chapter 5- We consider boiuidary-layer flows which are 
nearly similar, i.e, such that they are nearly described by 
solutions with one independent variable 7]!''. The limiting similar 
solutions corresponding to the prescribed pressure gradient 
parameter and second-order parameters are assimaed to be known 
for cr-«>. Corrections are sought in terms of the inverse 
powers of cr so as to extend the validity range of the solutions 
in the upstream direction. Alternately, such a problem may be 
posed as an initial value problem: an initial profile deviating 
slightly from the similar profile is prescribed at a given 
streamwise station, o\ and we wish to know how the similar 
solution is approached as cT-oo. Thus, we wish to confirm, in 
a mathematical sense, oirr intuitive expectation that the 
deviations of the initial profile from similarity will decay 
with increasing streamwise distance. 
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As already mentioned, various numerical and momentum 
integral methods are available to deal with nonsimilar boundary- 
layer flows and may be used to obtain boundary-layer developement 
downstream of a specified initial profile. However, we prefer 
an anal3rfcical method to formulate the problem outlined in the 
preceding paragraph. This method avoids the complexity of the 
numerical methods and is more accurate than the momentum inte- 
gral methods. According to the method employed here, the 
Variables are expanded in terms of a complete set of eigenfunc- 
tions multiplied by imspecified coefficients: these coefficients 
are to be determined from the initial condition. Although exact 
in principal, the method yields approximate solutions due to 
truncation of the series after a few terras. Such an approach 
has been used by Libby and Fox (1965) and Chen and Libby (1968) 
for the first-order boundary-layer flows. Further, Libby (l970) 
has given a review of the method (subsequently referred as Libby 
method) and its applications to a variety of problems. 

Following Libby method, the second-order nonsimilar 
boundary-layer solutions for the momentum transfer problem are 
constructed for large values of the streamwise coordinate, cr' . 

An initial profile (upto second-order) deviating slightly from 
the limiting similar profile ((r'-*<») is assumed to be known at 
a given station The solutions are linearized about the ' 

limiting similar solution and successive nonsimilar correction 
terms are obtained. The analysis leads to eigenvalue problem, 
the governing equations of which have a form similar to the 
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well~lcnovm Sturn-Liouville eq.uation and, therefore, it is 
possible to predict, in certain cases, that the eigenvalues 
are real and discrete; in particular, for Falkner-Skan flov/s, 
the eigenvalues are known to be real and discrete for the upper 
branch solutions in the usual range of p, i.e. 1988. . .<^<2 
(see Chen and Libby 1968). The solutions for the second-order 
effects associated with the longitudinal curvature, transverse 
curvature, displacement and external vorticity are obtained for 
a variety of parameters and are presented graphically. 

The results of the above analysis provide a ftirther 
insight into the singularities of the second-order solutions. 

The results of the first problem, described earlier, fix the 
location of the isolated singularities whereas those of the 
last problem show two distinct domains of stable and unstable 
solutions separated by the first critical value (corresponding 
to the singularities for various second-order effects mentioned 
in the first problem) of the parameter. Thus, it is shown that 
if the second-order parameters assume a value less than the 
critical value, any departure from strict similarity grows in 
the downstream direction and as such the flow is, in some sense, 
unstable , 

The application of the above procedure is illustrated 
by considering a uniform flow past a blunted wedge. The boimdary 
layer on a blunted wedge corresponds to a two-dimensional 
stagnation point flow (p = l) at the start (cr=:o)and then 
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accelerates to approach a constant (depending on the wedge 

angle) as Thus the initial and the final solutions are 

similar while at the intermediate stations , the flow is non- 
similar. A Gortler series solution, starting from the stagna- 
tion point, is already known from the second problem studied 
here. This series solution, extended sufficiently downstream 
by Eulerization , is used to prescribe the initial condition 
for the present problem for constructing the asymptotic solutions 
to give the approach to the far downstream (cr-o®) similar solu- 
tion. In addition to the eigenfunctions, we also include teimis 
arising due to variation of the second-order parameters with o ' , 
The multiplicative coefficient of the first eigenfunction is 
determined for each of the second-order effects, by joining 
the present series (cr* -*■«>) with the Bulerized G-5rtler series. 

The resulting solution provides in a certain sense a proof of 
convergence for the Gortler series solution. Further, the 
resulting asymptotic solution for large cT , in conjunction with 
the Eulerized G-ortler series may be used to assess the accuracy 
of the corresponding results obtained by certain numerical or 
approximate methods, where the error bounds are unknown or 
uncertain. 

It seems appropriate here to make a brief mention of the 
second-order theories other than that of Van Ifyke (1962). 
lenard (1961), Maslen (1963) and Honda and Kiyokawa (1969) have 
independently given the* second-order boundary-layer theory, the 
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first two authors obtaining solutions valid near the stagnation 
point while the latter authors formulated the problem in the 
stream function coordinates. More recently, an alternate 
approach to incorporate the second-order effects due to longi- 
tudinal curvature, displacement and transverse curvature in the 
first-order boundary-layer equations itself has been given by 
Van Dyke (1969), Davis, Werle and Wornom (1970) and Davis, 
Whitehead and Wornom (1971)* Werle and Wornom (1972) have applied 
this theory to study the longitudinal curvature and displacement 
effects for a circular cylinder in the incompressible flow. The 
single set of equations are suitable for numerical methods such 
as finite difference and finite element. However, since the 
second-order effects have not been linearly separated out, one 
has to solve the equations separately for each value of the 
parameter. Thus, for the present work where we study the contri- 
bution of the various second-order parameters to the individual 
second-order effects. Van Dyke’s (l962) formulation is preferred. 



CHAPIER 2 


FUNHAMENTAI. EQUATIONS AND RELATIONS 

2.1 Introduction : 

We describe here a procedure for obtaining the second- 
order incompressible bounda 2 ?y-layer equations, containing terms 
due to longitudinal curvature, transverse curvature, displacement 
speed, external vorticity and stagnation enthalpy gradient. Van 
I^ke (1962) has p 2 X)vided one such formulation for the momentum 
and heat transfer in constant density flow of a perfect gas. 
However, for liquids, the thermal problem of Van Dyke needs 
some modifications (as explained in Section 2.3). The main aim 
of this chapter is to present such a modified analysis for 
liquids. 

• The governing equations for the second-order boundary- 
layer problem are obtained from the Navier-Stokes equations, 
using the method of matched asymptotic expansions, with inverse 
powers of the square root of Reynolds number as the 

expansion parameter. This method divides the flow in two regions 
which are described by two different limit processes: an outer 
limit (defined as R-*<» with physical normal coordinate n fixed) 
valid everywhere except in a narrow region of order near 

the vmll and an inner limit (defined as R-«> with stretched 
physical normal coordinate nR*^'^^ fixed) valid in the region of 
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order near the wall# The repeated application of these 

limit processes yields outer and inner expansions which are 
matched in an overlapping region to obtain the missing boundary 
conditions. 

The coordinate system employed in the present work is 
described in Section 2.2. The Navier-Stokes equations are given 
in Section 2.3. A unified analysis for liquids and gases is 
presented in Section 2.4 and the final governing equations of 
the second-order boundary-layer problem are listed separately 
for liquids and gases in Section 2.5. These equations are the 
starting point for subsequent chapters and are studied for 
variety of flow situations. In the last Section 2.6, we 
also include a note on the displacement thickness which is 
needed in the subsequent chapters. 

2.2 Coordinate System ; 

The orthogonal coordinate system (s, n) used for the 
present problem is shown in Figure 2.1. Here s is the distance 
along the surface measured from the stagnation point or the 
leading edge and n is the distance normal to the surface. The 
velocity components in s and n directions are, respectively u 
and V. The longitudinal curvature of the body (reckoned positive 
for a convex surface) is denoted by k(s)* In axisymmetric flow, 
0(s) is the angle that the tangent to the meridian cuirve makes 
with the axis of the body and r(s) denotes the normal distance 
from the axis. The three quantities k, 0 and s are connected 
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tbroijgh the following relations: 

sine = dr/ds, cos 0 = -K"^ d^r/ds^. 

The length element dl in space is given hy 

dl^ = (1+Kn)^ ds^ Hh dn^ + (r+ncose)^^ d0^ 

where for the axisymmetric flow, j is unity and 0 is the azi- 
muthal angle while for the plane flow, j is zero and 0 is the 
linear distance normal to the plane of flow, 

2.3 Equations of Motion : 

The Navier-Stokes equations for a steady, plane or 

axisymmetric flow of an incompressible fluid in nondimens ional 

form (distances ^(a,n) are nondimens ional ized with respect to 

the characteristic length 1 of the body, velocity IJ(u, v) by 

reference velocity U., pressure by Pu^ and temperature T by 

c o 

div 2*0, (2.1) 

U.gradU +grad P = -R~^ curl curl U, (2.2) 

U.gradT - (PrR)'"V^ = R"''^ grad U.def U - eU. grad P (2.3) 

Here we have defined the Reynolds number R=pU L/j[iand the 
Prandtl number Pr =jaCp/k where p is the coefficient of viscosity, 
pthe density, k the thermal conductivity and Cp the specific 
heat. 

The quantity e, defined as (91ogp/91ogT)p, is the thermal 
expansion parameter. It may be noted that for a perfect gas its 
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Value is minus unity (e=:-l) while for liquid., it is negligibly 
small (e«0). The formulation of the thermal problem of Van 
I)ylse (1962) corresponds to e = ~1. In this case, the last term 
of the energy equation (2.3) can be eliminated with the help of 
the momentum equation (2,2) to obtain the equation 

D.grad (T+u2/2)_r-'' v^{Pr-% +-a^/z) = R" grad U).( grad U)*( 
used by Van Dyke (1962). However in the case of liquids, by putt- 
ing e=0 in the equation (2,5) one obtains the following equation: 

lJ.gradT-R“V^(Pr’"‘’T) R-’^ gradU-def U. (2.5) 

Since equation (2.5) does not follow from equation (2,4) as a 
special case, we have to reformulate higher order problem for 
equation (2.5), valid for liquids. We choose to formulate the 
problem for arbitrary value of e such that the governing equations 
for two special cases of gas and liquid will follow by choosing, 
respectively e=-1 and 0. 

The boundary conditions for the equations (2, 1)-(2.3) are 

U = 0 at n = 0, (2.6) 

T = T^(s) or • 6T/5n =0 at n - 0 (2,7) 

where (2,6) is the well-known no slip condition for impermeable 
surface. In (2.7), only one of the two conditions on the surface 
is required. The first prescribes the wall temperature T^(s) 
and the second implies an adiabatic wall (no heat transfer at 
the wall). As regards the upstream condition, the oncoming stream 
may be nonuniform (but prescribed, independent of Reynolds num- 
ber) and have gradients of velocity (vorticity) and temperature. 
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H.QJ^^iPUlatioD of the Se cond-Order Bo’undarv-.Layer Problem s 

We now obiiain bbe as3rnipto'fcic behavioiir of Navier— Stokes 
ec^uations for large values of Reynolds number using the techni- 
que of the matched asymptotic expansions. Application of the 
outer and the inner limits to the Ravier— Stokes equations yields 
two sets of equations which are matched in the overlapping region 
by using the well-known matching principle. 

2*4.1 Outer Expansion ; 

We write the outer expansion of the variables in equation 
(2.1)-(2.3) by considering the outer (Euler) limit defined as 
R -.00 with X fixed, The appropriate sequence of gauge functions 
for unseparated flows past analytic bodies consists of the 
inverse power of square root of Reynolds number (Van Dyke 1962). 
The outer expansions may be written as 

tJ(x; R) = ^ ^ ^ ^ ^ 

P(x;R) = P (x) + R~''/2p (x) + ..., 

Wo ^2.8) 

T(x; R) = T^(x) + + .. . , 

>i^(x; R) =^^(x) + + • • • * 

Here^ is the stream function, defined by. 

(r+ncose)^U= a>$/an, (l+Kn)(r+n cos e)^ V = -a'^as. 

On substituting the expansions (2.8) into the governing equa- 
tions (2.1)-(2.3) and collecting the coefficients of the like 
powers of gives the problems for the successive orders. 

The first order equations are 
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The first order equations are 


div U. = 0, 

(2.9a) 

U^.gradU^ + gradP^ ® 0, 

(2.9b) 

H • grad .( T ^ "feP ^ ) = 0 . 

(2.9c) 

These are the well-known inviscid Euler equations. The 
tion of these equations along a stream line leads to 

integra- 

+5^/2 = B^(^|/^), 

(2,10a) 

T^- eTj2/2= H^(^|^,). 

(2.10b) 

The Bernoulli function B^(>|^) and temperature function 
are to be evaluated from the upstream conditions. Equations 
(2.10a,b) show that the total pressuce and temperature are 


conserved along the first-order stream line = constant) 
defined as 

(r+ncose)^U^ = 

(1 + Kn) (r +n cos 0)^ = -S\J<|/ds. 

The second^order equations are 

^^2 =0» (2.11a) 

U^.gradlJg +U2.grad TJ^ ^gradP^^ 0, (2. 11h) 

U^.grad (T 2 + eP 2 ) +IJ 2 .grad (T^ + eP^) = 0. (2.11c) 

These equations are the perturbed Euler equations which repre- 
sent the displacement of the outer flow due to the first-order 
boundary-layer. The first integral of equations (2.11b) and 
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(2.11c) gives 

I’ 2 +Si*U 2 =xI^2®1^'^1^ +B 2 (>i<,), (2.12a) 

^2 — ’^2 ~ ^2^1 ^^1 ^ ^ ^ * (2.1 2b) 

where \|/2 is the second-order stream fimction for the outer 
flow and is defined by 

(r +n cos e)^ 02 = 

(1 +Kn)(r +n cos 0 )^ Y 2 = -dvj^/0s. 

As the oncoming stream is prescribed, independent of Reynolds 
number , we have ^2 “ ^2 ~ 

It is well-known that the outer equations do not satisfy 
the no slip boundary condition due to the loss of the highest 
order derivative of the Navier-Stokes equations at large 
Reynolds number. These equations are, therefore, not valid 
within a distance of order near the wall. It is clear 

that a different inner limiting process valid near the wall, 
in a region of order is needed, 

2.4.2 Inner Expansion ; 

To study the behavior of the flow near the wall in the 
region of order we introduce the stretched inner (Prandtl) 
variable 

1 = R'^/^n 

and study the limit R with N fixed. The inner expansions 
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are of the following form: 

u(s,n;R) = u^(s, N) + R"*^/^U 2 (s, N) + . . , , 

'v(s,n;R) = R”^/^ (s, N) + R“^ ^2 (s, N) + . . . , 

P(s, n; R) = p^(s, N) + R“'^/^P 2 (s, N) + . . . , (2.13). 

T(s,n;R) = t^(s, R) + R"''*/^ t 2 (s, N) -f- . . . , 

>i^(s,n;R) = (s, R) + R-'* ^ 2 ^®» . 

Substituting the above expansions in the governing equations 
and collecting the coefficients of the like powers of 
we obtain the problems of successive orders, The terms 
independent of give the following first-order equations: 


(r^u^)g + (r3v^)jj = 0, (2,14a) 

^l^ls ’*''^1’^1R“ '^1RR ”^^18 (2.14b) 

P^U = 0, (2.14c) 

with the corresponding boundary conditions 

u.j(s,0) = 0 , v^(s,0) = 0, (2.15a,b) 

t.j(s, 0) = T^(s) or t.jjj=0. (2.15c) 


Here subscripts s and R denote the partial differentiation. 

These are the well-known classical equations of Prandtl boundary 
layer. Again, the terms of order R~^^^ give the second-order 
boundary-layer equations as 
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[r^Cug + Nuy j cos 0/4?)] g + 3^^ ["^2 + 3 e/r)Fv^ = 0 ( 2. 16a) 

^2^s +^2^1N +’"l^2s - ^2M + ^28 

= K[(Nu^jjjr)jj -- v^ (lJu^)|j Hh j cos9/ru^jj5 (2*161)) 

^2N ^ > (2# 16c) 

Vls +''2*1]j "I'^as +^1^211 - *2NH 

= (K + 3 cos e/r) - 2Ku^u^j^ h- 2u^N^21T 

+ KNu^t^g- eMu^ p^g + ‘^'^1^28 ‘^^2^18 ‘^’^2^111^ ^ 

with the hoiandary conditions 

U2(s,0) = 0, V2(s,0) = 0 (2.17a, h) 

t2(s,0) = 0 or t2jj(s,0) = 0. (2.17c) 

2.4*3 Matching ; 

The outer expansion violates the "boundary conditions on 
the velocity and temperature at the wall while the inner expan- 
sion does not satisfy the boundary condition far away from the 
wall. Therefore, the validity of the inner expansion is res- 
tricted to a region of order near the wall whereas the 

outer expansion is valid only outside the region of 0(R''^'^^) . 

Thus neither expansion has sufficient boundary conditions and 
the missing boundary conditions have to be determined to make 
the problem determinate. This is done by matching the inner 
and the outer expansions in the overlapping region by the 
well-known matching principle (Kaplun 1957, Van Dyke 1964b) 
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written as follows: 

Inner limit of the (outer expansion) 

= Outer lim.it of the (inner expansion). (2.18) 

This matching principle is applied to all the variables and 
corresponding matching conditions are obtained. Since the 
procedure is identical for all the variables, we illustrate 
it for one typical variable say 0(s,n;R). The outer expansion 
is given by 

0(s,n;a) =<|^^(s, n) + n) + • •• . (2.19) 

To find its inner limit, we write the equation (2.19) in terms 
of the inner variables (s, R) and take the limit R with inner 
variables fixed. Thus, 

+ [<l>3(s,0) + N(f)2jj(s, 0) + 0)]+ ... _ (2.20 

The outer limit of the inner expansion is found by rewriting the 
inner expansion in terms of the outer variables (s, n) and taking 
the limit R with (s,n) fixed. This gives, 

0(s, nR”*/^; R-~) = 0.,(s, I--***) + (s, R-*co) 

+ R~V3(s, R . . (2.21) 

Using the matching principle (2,18) to compare equations (2.20) 
and (2*21), we obtain the following matching conditions: 

0.(8, R) =(|).j(s, 0) , 1 (2.22a) 

' ' -.oo 

02(s, R)'-R<|).j^(s, 0) +<|)2(s, 0). J 
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Using equation (2.22a) along with relations (2.10a,h), the 
matching conditions for the first-order boundary-layer problem 
may be written as follows: 

V^(s,0) = 0, (2.23a) 

u^(s,N) = U^(s,0), (2.23b) 

p^(s,N) = P^(s,0) = B^(0)-U^(s,0)/2, p-=o (2.23c) 

t^(s, N) = T^(s, 0) = H^(0) +eu2(s, 0)/2. ^ (2.23d) 

Equation (2.23a) represents the boundary condition at the wall 

for the first-order outer problem; the condition U^(s,o) = o 
has to be given- up. 


The matching conditions for the second— order boundary- 
layer problem are obtained by using equation (2.22b) alongwith 
(2.13). Thus, 


U2(s,N) v- U2(s, 0) + NU^^(s, O) = U2(s, 0) 


PgCsj N) 

t2(s,N) 


^^(s,N) 


+ N|;r3B:|(0) - K0^(s,O)] 

V. P2(s,0) + NP^^(s,0) = B’(0)Ni^2(®»0) 

- U^(s, 0) U 2 (s, 0) + KlU^(s, 0) 

V- T2(s, 0) + N T^^(s,0) = H'(0)^i'2^s,0) 
+ eU^(s, 0) U 2 (s, 0) +NU^(s, 0) 

[r^H^(O) - er^ B*(o) +eiaj^(s, 0)] 


) 


U -► oo 


(2.24a) 


( 2. 24b) 


(2.24c) 


(2.24d> 
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In the above equations (2. 24a-c), the extreme right hand side 
expressions are obtained with the help of equations (2.10) and 
(2.12). The matching condition (2.24d) can be rewritten in the 
following forms: 

>i^.(s,0)= lim (1/^. -Nfc), (2.25) 

O) = lim (v.| — NVjjj*)# (2«26) 

JJ oo 

This is an important matching condition for the second-order 
outer flow and represents the effect of displacement of the 
inviscid (outer) flow by the boundary-layer (inner) flow. The 
quantities on the right hand side of the expressions ® 
and ¥2(3,0) can be evaluated from the first-order boundary- 
layer solution, A physical interpretation of these equations 
is that the effect of the classical boundary-layer upon the 
outer flow is equivalent to surface distribution of sources, 

2.5 S ummary of the Governing Equations ; 

2.5,1 Equations for the outer flow : 

Eirst Order (Euler Elow) : The governing equation of the first- 
order outer stream f unct ion , obtained from equations (2,9) 
and ( 2 . 10 ) is 

= r^(l +Kn) B^(l'^) (2.27a) 

with the corresponding matching and boundary conditions given by 

’^^(s,0) = 0, ="^100 far upstream (2.27b, c) 

where and’^^^ are known from the upstream conditions. 
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The integration of equation (2»11a) subjected to the boundary 
conditions (2.11b,c) gives the velocity components U^(s,n) 
and y^(s, n) and, in particular, U^(s, o) needed in the matching 
conditions for the first-order boundary-layer problem, Further, 
using and V^, pressure distribution can be determined 
by the Bernoulli equation (2.10a). 

Second Order (Displacement Flow): The governing equation of 
the second-order outer stream function are obtained from 
eqi^ations (2.11) and (2.12) and can be written as: 

(r“3 (1 + Kn)-1) (r"*^ ( j + Kn) ^ = 

r3 (1 + En) b" (v|/^)^g , (2.28a) 

with the corresponding boundary and matching conditions as; 

\p2(s,0)= lim(^-Bfe), 

^ -►oo * *■*' 

= 0 far upstream. (2,28b, c) 

where is obtained from the first-order boundary-layer solu- 
tion. The solution of the problem (2,28) yields U 2 (s,n) and 
whereby we can evaluate IJg(s, O) required for the 
matching conditions of the second-order boundary-layer 
problem. Using Ug and Yg alongwith equation (2.12), we can 
evaluate the pressure distribution Pg. 

2,5.2 Boundary-Layer Equat ions ; 

We now write down the governing equations of the first 
and second-order boundary-layer problem for both liquids and 
gases. Since the continuity and momentum transfer (velocity 
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boimdary-.layer) eq.'uations are same for both the cases 
(liquids and gases), we write them once only. The energy 
(thermal boundary-layer) equation for the liquids and gases 
is obtained by putting respectively e=0 and 1 in the 
general form of the energy equation, boundary conditions 
and matching conditions obtained earlier in § 2.4- and are 
given separately for these two cases. 



Birst-order-' 


(r^u^)g +(r?v^)jj = 0, 

( 2. 29a) 

s ■^'^1^1n “ 0) U^g(s, 0) , 

( 2. 29b) 

u^(s, 0) = v^(s, 0) = 0, 

(2,29c) 

u^(s,]lir) !=ll^(s,0) as , 

( 2. 29d) 

Second-order: 


+Hu^ 3 cos e/r)] g + [vg + (K + j cos e/r)Nv^] jj 

= 0, (2.30a) 


UlUja+HaUis+ViUjjj+VgUuj-UajH = 

- m^{3, 0) 0) + - [Kinj2 (g, 0) 

oo 

+ K / (U^(s, 0) - u^) <aN]g +u^^ 3 cos 0/r 

- (v^-Nv^)] +[U^(s,0)U2(s,0)]g , (2.50b) 

U 2 (s, 0) = 0 =V 2 (s, 0) , (2.30c) 

U 2 (s,N)-' - (s, 0) +r3HB^(0) + 0) as . (2.30d) 
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M3. m y ^e.quat ion for liauids 

W II 111' 1 1 I III B HI ni"lll■^llll Br'i«iiiu.i 


Pirst-order: 


"in =0' 

(2.31a) 

t^(s,0) = T^(s) or t-ijj = 0, 

(2.31b) 

t^(s, N) =.H^(0) , as N-oo . 

(2.31c) 

Second~order: 


Uit23+U2ti3 + 2u,1,U2N 


= K[Nu^t^^+Pr-''t^jj- 2u.,u^^] + Pr-'' j oos 

0/r, (2.32a) 

t 2 (s, 0) = 0 or t 2 jj(s, O) = 0, 

(2.32b) 

tgCs,!!) - [vl/^Cs, 0) + r^NU.^(s, 0)] h'( 0) asN-oo. 

(2.32c) 

Energy equation for gases fe=-1'J: 


Pirst-order: 



(2.33a) 

b^(s, 0) = t^(s) 

(2.33b) 

h.^(s, N) = E^iO) as N-oo . 

(2.33c) 

Second-order: 


^1^2s ■’■^2^1s ■^'^1^2N ■'*’^2^1N ^2NW 


= K[Nu^h.jg +Pr“’' h^jj] + Pr""* h.^^^ i cos e/r, 

(2.34a) 

hgCs, O) = 0, 

(2.34b) 

hgCs, IT) - [Vi/2(s, 0) + r^ NU.j(s, 0)] H^(0) asN-“. 

(2.34c) 


It may be noted that the above form of the energy eq.ua- 
tions for gases is preferred in order to conform it with that 
of Van Dyke (1962). It is obtained by substituting e=~1 in the 
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energy equations (2. 14-d) and (2,16d), which is then manipulated 

with the help of the momentum equations to obtain the above 
form. 

2.6 Displacement Thickness : 

Before closing this chapter we define and obtain expres- 
sions for 6 and - the first and second-order displacement 
thickness. The displacement thickness 6 is defined by the relation; 

^OO ^ ^OO 

f Ur^dn = / ur^dn ( 2 . 35 ) 

where n^^ is the value of n arbitrarily chosen to be beyond the 

outer limit of the inner layer. The displacement thickness may 
be expanded as; 

5(a; R-1/2) = 11 - 1/25 ( 3 ) +R-1 A (s) . . . . (2.36) 

No?^, we letR-*ooas n^^ remains fixed and rewrite equation ( 2 . 35 ) 
in the following form: 

6 .. 

/ U(s,R-‘'/2u) r3 dN 

0 

oo 

= f [u{s, R“^/^ H) - u(s, IT)] r^ dN ( 2 . 37 ) 

The expansions for u(s,F) and XJ(s,R ^^%T) are given by respec- 
tively, equations (2.8) and (2r13). Also we have to expand r^ as 

r3 = + R-1/a H 3 cos 0 (2.38) 

Purther, the expression for D(s, R'^^H) is expanded in a laylor 
series about R = 0 so that we have 

u(s,R-''/%).-n^(s,o) + R-V2 [;]ju^^(g^o) ^ ^ 

Subs'fc iiiu'feing "the above mentioned expansions in ec[uaiion ( 2 # 57 ) 
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and integrating the left hand side, we obtain: 

U^(s,0)r3 6 + [u^(s,0 )r^A + ^ (s, 0) 

+ 6 U^Cs, 0) r^ +-y U^(s, 0 ) j cose] = / [lU^(s, 0) - U^(s, N)}r^ 

+ N[U^(s, 0) - u^(s, N)] 3 cose]] dN (2*40) 

Now collection of the coefficients of the like powers of R-Va 


yields the following expressions for the first and the second- 
order displacement thickness: 


6 =/ [1- 


0 


A = - 


u.(s,N) 

07 ( 3 , 11 ) 3 “ > 


(2.41) 


1 o2 loose r u.(s,N) , 

2 ^i(s,0)“ - r. 2 ’ ^ I'l J 

Miller;. A. - ^U,,^(s,0), 

uY(s,o) !■» irr£o) (2.42) 



CHAPTER 3 


SEIF- SIMILAR SOLUTIONS OF THE SECONJP-ORDER BOUNLARY-LAIER 
EQUATIONS INCLUDING DISSIPATION TERIiS 

3. 1 Introduct ioa ; 

We now study the self-similar solutions of the second- 
order houndary-layer equations including the dissipation terms, 
viz. the equations (2. 29)--( 2.34) given in ^2.5.2 of Chapter 2. 
In these equations, the terms involving K arise due to longi- 
tudinal curvature, 3 cose/r due to transverse cuirvature, B^(o) 
due to external vorticity, H^( 0 ) due to stagnation enthalpy 
gradient and the remaining terms due to displacement speed. 
These effects are of order as compared to the classical 

houndary-layer theory and become significant when R is only 
moderately large,. 

The above mentioned second -order problems have been 
studied by various authors and an excellent review is given by 
Tan Dyke (1969). Recently, Afzal (1969, 1972) a,nd Werle and 
Ifevis ( 1970 ) have proposed a more general class of similarity 
as compared to the earlier works and have shown that the des- 
cription of the second-order problem needs three principal 
functions, namely those due to longitudinal curvature, trans- 
verse curvature and displacement speed, Af^-al £ind 0berai(1972) 
have studied the momentum and heat-transfer in the second-order 
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■boundary-layer flows without the viscous dissipation terms 
while Werle and Davis (1970) have studied only the momentum 
transfer. Afzal and Oberai (1972) have shown that for a given 
pressure gradient parameter, the solutions to the second-order 
equations show a sequence of infinite discontinuities for cer- 
tain critical values of the second-order parameters while 
Werle and Davis (1970) have reported only one such discontinuity. 
However, these authors have not given a detailed account of the 
o-ccurrence of such discontinuities. We investigate here (^3.3) 
the general pattera of these infinite discontinuities and suggest 
their relevance to the design of body shapes so as to achieve 
some desired aerodynamic properties. 

The main aim of this chapter is to study the general 
self-similar structure of the second-order boundary-layer equa- 
tions when the viscous dissipation terms are included. The con- 
tribution due to the viscous dissipation may be significant 
whenever the product of the Prandtl number and Eckert number is 
of the order of unity. We may expect one such situo,tion for 
liquids of moderately large Prandtl number. We consider here 
steady, laminar flow v/ith constant density, past a two-dimen- 
sional or axisymmetric body. The study includes both types of 
the thermal boundary conditions at the wall: ( i) specified 
wall temperature and (ii) an insulated wall (no heat transfer 
at the wall). Eor the latter case, the classical concept of 
the recovery factor is extended to the second-order flows. The 
terms in the heat transfer and recovery factor due to 
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longitudinal curvature, transverse curvature, displacement 
speed, external vorticity and stagnation enthalpy gradient are 
calculated for various values of the parameters arising in the 
problem and the results are displayed graphically. 

We first describe in §3*2, the momentum transfer problem 
and then in §3*3, discuss the occurrence of the infinite discon- 
tinuities. The subject matter described thereafter is sub- 
divided into two parts. In part I (§§3.^-3*7)j we study the 
theimal problem for liquids while part II (663.8-3.11) is 
concerned with the corresponding thermal problem for the cons- 
tant density flow of gases. 

3.2 Momentum Transfer Problem : 

The momentum transfer problem is governed by equations 
(2.29) and (2.30) given in ^2.5.2. In viev/ of the reasons 
given in Chapter 1, we choose the Gortler variables to study 
the boundary-layer equations for nonsimilar flow situations. 

The Gortler variables are defined as 

s 

0 -= / U^(s,0) r^3^s , 7? = r3(2cr)“'’u^(s,0)H. (3.1) 

G 

The equation of continuity (2.29a) for the first-order 
boundary-layer is identically satisfied by introducing the first 
order stream function ^ such that 

. (3.2) 

It is found convenient to write the stream function in the 
following form: 

= V(2o') 


(3.3) 
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Introducing the above expression (3»3) in the first-order 
momentum transfer equation (2.29), we get, 

f ff »' +y\u^(l-f '2) - 20" (f 'f^ -f ”f^) = 0, (3.4a) 

f(cr, 0) + 2cr f^(r, 0) =0, (3.4b) 

f'(cr, 0) =0 , f(cr,oo) =1. (3.4c,d) 

Here prime ( ') denotes the differentiation with respect toT^. 

The quantityAy^ is the principal velocity function defined by 

= t2o'/U^(s, 0)] [dU^ (s, 0)/da'] . (3.5) 

The function A is important parameter for the boundary- 
layer problem. It has relevance for nonsimilar as well as similar 
flows. Por the latter case, it remains constant in the streamwise 
direction and can be identified as the familiar Palkner-Skan 
pressure gradient parameter denoted by p. However, for nonsimilar 
flows, A varies in the streamwise direction (see, e.g. the 
blunted wedge problem considered in Chapter 4). In general, Ay^ 
can be determined, if the first-order outer velocity, U^(s,0), 
is known either from the solution of the outer flow or by curve 
fitting to the experimental data. 

The continuity equation (2,30a) for the second-order flows 
is satisfied by defining the second-order stream function Xl /'2 Q-® 

[u^ + dNu^ cos 0/r ] = r^ [v^ + (K j cos 0/r) Hv^] = - (5, 6 

If we express the second-order stream function in the following 
form 


^2 = f(2cr)T'(cr,7/) 


(3.7) 



35 


the second-order momentum transfer problem ( 2 . 30 ) reduces to 
+ 2 Au^ f 'I" +f' *P - 2cr (f^, p ' _f 

= l)(t + + 2 cr f 4cr/f^f»d7;; 

+ (Abi + +a + 2<Ta^)]/(l +J\^^ ) ^ 

+ +f "') 4-f " + ff ' -Aj 3^ f ' 2 + 2crf .f 'j 


B^( a + 2cra^) _ B^( 2Ag^ ’^Ab^) > 

P(cT-, 0 ) + 2 o'P^((r, 0) =0, P‘((r,0) =0 , 

P’(o,~) =:7;(^b^+B^ +B^ . 


( 3 . 8 a) 

(3.8b,e) 

(3.8d) 


In the above equations, we have defined the following quantities 
for convenience of writing: 

’ Bp = f(20') K/r^ U^(s, 0) , (3.9a,b) 


= f(2cr) jcose/r23u^(3,0) , B^ = U2(s, 0)/U^(s, O) , (3.9c,d) 
B^-= f(2cr)B^(0)/u2(s,0) . 

Purther, the quantities A b^ » A and Ab^ are respectively the 
principal functions corresponding to longitudinal curvature, 
transverse curvature, and displacement speed and are defined by 
Ab^ = ( 2 cr/Bp(dB./dcr) , i = l,t,d (5.IO) 

It is noted that the external vorticity does not give rise to 
a new principal function since 


- (2cr/B^) (dBydcr) = 1-2^ . 


( 3 . 10 a) 


Again, theAB’s are constant for snnilar flows while for non- 
slmilar flows, they vary in the streamwise direction (see, e.g. 
the bliinted wedge problem in Chapter 4). 
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In the present chapter, since we are interested in 
studying the general self-similar solutions, the modified 
stream functions tier ,V) and 1(0", 7)) are assumed to he of the 
following form? 

f(<r,7)) - (3.11a) 

l'(cr,?7) = Z i = l,t,d,v (3.11h) 

i ^ 

In writing the eq.uation (3.11b), v/e have exploited the linearity 
of the second-order problem and decomposed the contributions of 
various second-order effects which have a clear physical inter- 
pretation. 

Substitution of (3.11a) in (3.4) yields the following 

ordinary differential equation for the first-order boundary-layer: 

f »*'+ ff ” + 6(1 - f '2) =0, , , 

( 3 . 12 ) 

f(0) = 0 = f ’(0) , f'(“) = 1 

which is the well-known Balkner-Skan equation. Here the princi- 
pal velocity functionA^^ is constant and is denoted by-p. 

The second-order boundary-layer equations, obtained by 
substitution of (3.11b) in (3.8), reduce to the followiag system 
of ordinary differential equations; 

f «««(i) + fp>t(i)_ (2P + X.) f 'I"^^)+(1 +Ai)f 

f .. ,.s „ i=l,t,d,v 

(O) = 0 = P (O) , (Tj) = GjL as V 

(3.13) 

where ^i given in Table 3.1. Here again, the 

principal functions A 3^ ,A and are constant and are 

respectively denoted byA]^, A-^ andA^ . 
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Table 3.1i Values of A's, b's and C's in the second-order 
boundary-layer equations. 


Second-order 

effect 

i 

Ai 

L. 

1 

Gi 

Longitudinal 

1 

^1 

^i)(f" + ff') 


curvature 



+ (2p+A]_)(p77+a)]/(l + p) 

(3.14) 

Transverse 

curvature 

t 

^t 

-.■J7(2p + f •”) +f '’+ff ' 

7? 

(3.15) 

Displacement 

speed 

d 

^d 

- 2p -A<1 

-1 

(3.16) 

External 

vorticity 

V 

1-2P 

- a 

7? 

(3.17) 


Solutions to the momentum transfer equations (3.11) have 
been adequately described by Y/erle and Davis (1970) and Afzal 
and Oberai (1972). Afzal and Oberai (1972) have shown that for 
a given p, the solution of the second-order equations contain 
a large number of singularities for some critics,! values of 
the parameters /\_i»A-t Davis (1970) also 

pointed out the first, of these singularities. Since no attempt 
is made by these authors to elaborate on the observed infinite 
discontinuities, we investigate, the general pattern of the 
singularities and point out their mathematical as well as phy- 
sical significance in the next section. 
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5«3 Investigation of the Singular Behaviour of the Second-Order 
Solutions for Some Critical Values of the Second~Order 
Parameters ; 

First of all, we wish to establish that the singularities 
in the second-order solutions occur whenever the parameters A 
and A assume a value equal to the eigenvalue of the correspond- 
ing second-order homogeneous problem given by 

I*(M, A .)“ + - (2P+AJ f 'M’+ (1 +AJ f ”M *= 0, 

^ ^ (3 18) 

M(0) = M'(0) = = 0, 

The operator (3»18) is similar to the one encountered by Chen 
and Libby (1968) in connection with the asymptotic studies of 
the first-order boundary-layer theory. Chen and Libby (1968) 
have estimated the first twenty eigenvalues corresponding to 
Poo = 0,5 and the first eigenvalue corresponding to several other 
values of p» 

The contribution of the longitudinal curvature, transverse 

curvature and displacement speed to the skin friction for two 
typical values of p = 0.5 and 1,0 is shown as a function of the 
corresponding second-order parameter in Figures, 3.1, 3.2 and 
3.3. The eigenvalues of the associated homogeneous operator 
(3.18) are taken from Chen and Libby (1968) and are indicated 
by an arrow on theA^-axis in the above mentioned figures 
(in, n = 1,2, 3, 4 denote the first foixc eigenvalues for p == 0,5 
and D^, the first eigenvalue for p = 1.0), Figures 3, 1-3.3 show 
that for positive values of the parameters A 2^, A ^ andA^^, the 
corresponding: second-order contributions to 'the skin friction 
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ar.-well,beha.vo4.. However, for the negative values of 
the parameters, a large number of stagularitles occur; the 
location Of these singularities is seen to be near the values 
Of the parameters A, andA, such that it is egual to the 
eigenvalue of the operator (3.18) indicated by arrows. This 
behavwur is further supported by the close fonn solution for 
the displacement speed problem, which Is given by 


^ = KP- l)7)f' +f]/P 

forAj = p-2. Solution (3.I9) shows that for p = 0 , there is 
an infinite discontinuity atA^ = - 2 ; Ohen and Libby ( 1968 ) 
have shown that for p = o, the first eigenvalue of the operator 
(3.18) IS indeed - 2 . Thus, the general pattern of the singula- 
rities in the second-order solutions can be predicted: the 


location of the isolated singularities is given by the values 

of the second-order parameters A 1 , A -t andAj such that they are 

eaual to the eigenvalue of the corresponding homogeneous 
operator 5.18. 


I’urther, it is well-known (cf. e.g. Courant and Hilbert 
1953) that the unique solution to the nonhomogeneous boundary 
value problem exists only if the corresponding homogeneous 
operator with homogeneous boundary conditions possesses trivial 
solution. However, if the homogeneous operator has an eigen 
solution, then, the nonhomogeneous problem does not have a unique 
solution; the bounded nonunique solution exists only if the 
eigenfunctions are orthogonal to the nonhomogeneous terms of 
the given differential equation. In the present study, the 
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solutions at the eigenvalues hecome unbotinded (in other words, 
do not exist) and, therefore, suggest that the associated 
eigenvalues are not orthogonal to the nonhomogeneous terms of 
the differential equations. 

T/e suggest that the above observations may be used as a 
caution while designing the profiles of aerodynamic bodies. 

In order that the contributions due to second-order effects to 
the skin friction remain boimded, the profile should be chosen 
such that the values of the second-order parameters A A ^ andA^j^ 
are not in the neighbourhood of the above mentioned eigenvalues. 

It seems that, when the values of these parameters A » A^ andA^j^ 
approach a value equal to eigenvalues, the skin friction etc. 

become very large and as such the second-order boundary-layer 
theory of Van I)yke (1962a) breaks down. 

We now cite an analogy (from solid mechanics) with the 
linear oscillations of a spring-mass system under an applied 
periodic force. Although the analogy may seem rather far-fetched, 
nevertheless it serves well in understanding the significance of 
the singularities discussed above. For a spring -mass system, it 
is well-known that when the frequency of the applied force 
approaches the natural frequency of the system, the amplitude 
of oscillations increases leading to the phenomenon called 
resonance (infinite amplitude). In designing such systems, 
sufficient care is taken to separate the two frequencies so that 
the operating frequencies are kept distinct from the natural 
frequencies. Further, the linear theory breaks down in the 
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neighboiirhood of resonance and a nonlinear -theory is needed 
(of. e.g. Cole 1968, Narasimha 1968). Analogously we suggest 
that in the neighbourhood where the second-order parameters 

andA^ assume a value close to the eigenvalues of the 
associated homogeneous operator (3»18), a nonlinear boundary- 
layer theory seems in order. 

Part I: Heat Transfer in liquids 

3.4 Thermal Boundary-Laver Equations for Self-similar Solutions : 

The governing thermal boundary-layer equations (2.31) and 
(2.32) are given in §2.5 of Chapter 2. We shall again use the 
Gortler variables (o-',?)) defined by (3.1). Pirst, to obtain the 
general nonsimilar equations for the thermal boundary-layer , we 
write the first and second-order temperatures, t^ and tg, * 
in the following form; 

=[I„-H^(0)] g((r,7J) +H^(0) (5.20) 

tg = [I^-Hi(0)] G(<r,-^) (5.21) 

Introducing (3.20) in (2.31) > we obtain the first-order themal 
boundary-layer eqmtion in the following form: 

Pr“‘' g»' +fg' -Am H g - 2 cr (f'g -f g') = -Ef"2 , 

■^w”"l . (3.22) 

g((r, 0) = 1, g(cr', ~) = 0. 

Substitution of (3.21) in (2.32) yields the second-order thermal 
boundary layer equations. Thus, 

Pr-''G''+£GI..A5,^jj^f'G-2o-(f'G^-f^G') = I"g 

-Pg' - 2tr(P'g3,-P^g') - 2Bf"P" 

+ Bi[-Pr”n7?g') ' +Sf"('??f" +2f')] 

+ B^[-Pr'"1(7?g') • +Bf "(7)f "+2f‘)] . 

G((r, 0) = 0, &((r,T]f) = Bjj( 17 - a) as 


(3.23) 
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In the above equations, 
defined as 


A 




is the thermal principal function 


9[T^-H^(0)]/acr. (3.24) 

Ihe thermal principal function, Am u » will be constant for 

iw-^1 

self-similar flows whereas for nonsimilar flows, it varies along 
the streamwise direction. Por convenience of writing, we have 
introduced Eckert number, S, and a quantity Bjj defined as 

E (3,0)/[I^-H^(0)], (3.25) 

Bh = V(2<r) H^(0)/[T„-H^(0)] . (3.26) 

Since we are interested in studying the self-similar 
solutions, we assume that the fiinctions g and & are of the 
following form: 

g((r,77) = g.^C'^) + EggCT?) , (5.27) 

G(o-,7?) = S (7/)] , i=l,t,d,v,H (5.28) 

while Afp Ti is constant, denoted by A . In writing the above 
equations, we have exploited the linearity of the energy 
equations and separated the heat transfer due to the prescribed 
wall temperature [g.|(0) and G^(0)] and viscous dissipation 
[g^(0) and G2^0)j . Thus the heat transfer problem can be 
solved Independent of the Eckert number. Eurther, in writing 
(5.28), we have again subdivided the second-order problem into 
various physically meaningful problems; as was done for the 
momentum transfer case. Finally, it can be easily shown that 
A is the index of power law for wall temperature variation: 

T^ = T^ + Ccr”^/^, where C is a constant. 
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Subrtltutlon of ( 3 . 27 ) (j.gg) ^ 

yfelda the followtog ordtoary differential equations of the 
thermal boundary-layer problem: 

£irst-Qr(iP-r t 

e\' - f s\ A fg^ = 0, 

g^(0) = 1 , g.j(oo) = 0. 


Sg - fgg- 2pf 'gg = _f "2 ^ 

§2^®) = 0 = g2(°°). 

Second-Orr}p-P « 

Pr-^o';(i)tfa;Ci).(A.Apf'.(^).(i.A,).(i)g' 

-AF'(i) = M^, 

eS"\o) =0, (}(i)(7,) =a^ ,, 

Pr-^a;(i).fa'(i). (2p.Aa)f'a(i).(i,App(i)g: 

f / 4 N If /* • N 2 


(3.29) 


(3.30) 


(3.31) 


4^ho) = 0 = (oo). 


— JL' C ' 

- 2pF'(i)gp + 2f " F"(i) = If 

= a(i) 


(3.32) 


where Ai, M., and d^ for various second-order effects 

(denoted hy superscript i = l,t, d,y, H) are given in 
Table 3,2. 

It may be noted that even when the dissipation terms 
are included, the similarity is obtained when the Eckert num- 
ber, E, in constant; this leads tc A = 2p. This ccnditlcn 
implies that the wall temperature cannot be prescribed arbi- 
trarily; it Will depend on the value of the pressure gradient 
parameter for the given problem. 
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Table 3.2: The values ofA’s, d'e. M'e and Ifs In the second- 
order thermal boimdary-layer equations for 


liquids. 


Second-order 

effect 

i 

Ai d^ 

N. 

1 

Iidngitudinal 

curvature 

1 

Al -Pr-'’(7)g^') ' 0 

-Pr“ (TJg^) ' +f ”( 2 f 'j)^f " ) 

Transverse 

curvature 

t 

■^t 0 

''(i9g’) '+f"(2f'-?]f '•) 

Displacement 

speed 

d 

0 

0 

0 

External 

vorticity 

V 

1-2p 0 0 

0 

Temperature 
gradient in on- 

H 

^“A 0 7]-c< 

0 

coming stream 




Transfer and Recovery Eaetn-r« 


The total heat transfer rate at the wall is given by 

« = -«aT/an)„^, 

where k rs the thermal conductivity of the liquid. Substituting 

the expansion of T from (2.13) using relations (3.1), (3.20), 

(3.21), (3.27) and (3.28), the above expression for q may be 
written as 

q =-Pr-1r3n^(e,0)R-1/2(I^.Hp(2a-)-V2(,^^S-1/2^ 


= gj(o) -fEg^Co) 

q 2 = 2 Bi [g^ ( i) (0) + E a^( i) (o)] , 


1 ^ 


l»t,d,v,H (5.54) 


with 
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As mentioned earlier, the full similarity of heat transfer 
problem with dissipation exists only for A= 2p which corres- 
poiads to constant Eckert number. 

The temperature of an tooulated wall Tj, (also known as 

the recovery temperature) Is defined as the wall temperature 

for which the total heat transfer rate q is zero. The expression 

(3.3+) with q = 0 gives the following expression for the recovery 
temperature: 


= H^CO) “U^(s, 0 ) g 2 ( 0 )/g.j( 0 ) +u2(s, o) R“''/2 
f ®i[lg 2 ^ 0 ) G^(i)(0) -g'(0) a’(i) ( 0 )>/g’ 2 (o)] 

-R 0) G^’^^)(0)/g^(0) , i=l,t,d,v ( 3.35 

where Bg = f(2(r)H;(o)/u2 ( 3 , o) . A move convenient way to study 
the recovery temperature is to express it in nondmensional ' 
form, called the recovery factor and defined as 

rf = 2Ct^.-H^(0)]/u2(s,0) (5^5^. 


It may be noted that the recovery temperature and recovery factor 

are functions of Reynolds number R and also call for the following 
expansions as R->°o; 


with 


Tr = Tr.^ + R"''/^ 

= r^ +R-V2 j,2^. 
rg = S Bi r^i) 


. 


(3.57) 

i = 1 , t, d, V, 1 (3p58) 


Substituting (3.55) in (3.36) and using (3.38), wo obtain the 

following expressions for the first and second-order recovery 
factor: 
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(3.39) 


= -2g’(0)/g^'(0) 

4 2 [gjco) - 8 ^( 0 ) e’^^’(o)]/gj 2 ( 0 ), 

_(H)_ r^V,,, i = l,t,d,7 (3.40a) 

^ 2-2 r®-|( 0 )/g.(o)] 

(3.40b) 

Introducing the recovery temperature ( 3 . 35 ) into th^ 

(3.34) for the total heat transfer rate , and with A = 2B we 

0 tarn the following expression for q: 

q = -Pr~‘*r^ (P^-, n~ 1/2 p- 1/2 „ / , . 

, ’ R ^'l(s.O) [\-Tj 

[gl(0)+R-l/2 j j „r(l) . ■ 

i=l,..v 1 (3.41) 

^hls Shows. When the wall to.perature, equals the rooovery 
temperature, there is no heat transfer at the wall. 

Of he tT 

Of heat transfer, as 

«h =^(2<r) V2pr r 1/2/[;^S 0)( V^r)] 

= -8^(0) -R~ 1/2 2 3 £,'(!), , 

^11 lo;, 1 _ i,t,d,T. (3.42) 

msplaoewent 

The dlsplaoement speed problem for jointly self-slmUar 

°-e. i.e A, = 0 has the following olose form solution: 

= (7)f'+f)/2, 

r>(<i) / (3.43a) 

=7?g/2, 

r(d) ^ , (3.43b) 

^ (3.430) 

ustog the above solution, the seoond-order oontribution to 

reoovery factor can be expressed in terms of tbe first-order 
recovery factor by the following relation: 

4 ^’ = 2r.. 


(3.44) 
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5'. 7 Dlsousslort : 

The heat-transfer at the wall and recovery factor are 
obtained by Integrating equations (5. 29)-(3.32) along with those 
of momentum transfer, namely (3.12) and (3.13). The numerical 
integration was performed by using Runge-Kutta-G-ill method on 
IM-7044 Computer at Indian Institute of Technology, Kanpur. 

The step size used wasA*7= *05 and all the calculations were 
carried out in double precision to minimise the round off 
errors. The equations were integrated for full similarity 
with dissipation which corresponds to A = 2p. Contributions 
of longitudinal curvature, transverse curvature, displacement 
speed, external vorticity and temperature gradient (in on 
coming stream) to heat transfer, G^(0), and recovery factor, 
rg are displayed graphically. It may be noted that the heat 
transfer G^(0) [due to proscribed wall temperature] , ^^(0) 

[due to viscous dissipation] and recovery factor, are 
related to each other so that given any two of those quantities, 
the third can be determined. Therefore, we present results 
only for two of these quantities, namely heat transfer, G^(o), 
and recovery factor, r^, in the range -. 1988. .,< p< 2.0, 
0.5<Pr£5.0 and “2.0 <A]^> At >^d - shall first discuss 

below each of the second-order effect separately and then a 
comparative study of these effects will be given for a repre- 
sentative case. 

The second-order contribution due to longitudinal curva- 
ture to heat transfer G^(o) is shown in Figure 3*4« A large 
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ntan'ber of infinite discontinuities are observed forA^ < 0 v/hile 
£orAj_ > 0, the heat transfer is well behaved. Such types of 
discontinuities for < 0 are also reported by Afzal and 
Oberai (1972). An elaborate discussion of such discontinuities 
was given in §3.3 and, therefore, it is unnecessary to discuss 
results for negative values of A]_. PorA^ > 0, Pigure 3.4 shows 
that for given Ai and Pr, the heat transfer G^(o) decreases as 
P increases. Purther, the heat transfer g’(o) increases as A3_ 
(alternately Pr) increases for a given Pr (given Aj^) and p. 

Since the sign of the second-order heat transfer G^(o) is oppo- 
site to that of the first order heat transfer g‘(o), the positive 
longitudinal curvature (convex surface) will decrease the total 
heat transfer. The second-order recovery factor r^ due to 
longitudinal curvature is plotted in Figure 3.5 which shows 
that the recovery factor r^ decreases (in magnitude) as p increases. 
Purther, the recovery factor increases (in magnitude) as A i) 
(alternately Pr) increases for a given Pr (given A^) and p. As 
the first and second-order recovery factors are of opposite 
signs, the overall recovery factor decreases on a convex surface. 

The effect of the transverse curvature on the second-order 
heat transfer G^(0) and recovery factor are, respectively 
shown in Figures 3.6 and 3.7. Again we observe several infinite 
discontinuities forA^ < 0 which are of similar nature as for the 
case of the longitudinal curvature and, therefore, no additional 
comment is needed. PorA^ > 0, we observe from Figure 3.6 that 
the heat transfer G^(o) decreases (in magnitude) as p increases. 
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further, as A ^(alternately Pr) increases for a given Pr (givenA^) 
and p, the heat transfer G^(0) increases in magnitude. The sign 
of heat transfer G-^(O) being negative (which is same as that of 
the first-order heat transfer) , the total heat transfer will 
increase for positive transverse om?vature. Figure 3.7 shows 
that the recovery factor r^ decreases as p increases. For a 
given Pr and p, the recovery factor increases asA^ increases. 
Further, the recovery factor r^ is positive for Pr > 1 and nega- 
tive for Pr < 1 which leads to the conclusion that the total 
recovery factor for positive transverse curvature will increase 
for Pr > 1 and decrease for Pr < 1 . 

The solutions to the displacement speed problem for the 
heat transfer G^(0) and recovery factor are shora in Figures 3 »'$ 
and 3» 9 respectively. Here again we discuss the results for 
A(i>0 only for reasons mentioned earlier. Figure 3.8 shows that 

heat transfer G^( 0 ) increases (in magnitude) as 
(alternately Pr) increases for a given Pr (given A and p. 
Further, the sign of the heat transfer G^(o) being inegative, the 
total hea.t transfer will increase for positive displacement 
speed. From Figure 3. -9, we observe that the recovery factor 
r^ increases asAj^^ (alternately Pr) increases for a given Pr 
(given A^) and p. Exact solution for the case ofA^ = 0 given 
by (3.44) shows that the second-order recovery factor r 2 can be 
expressed in terms of the first-order recovery factor r.| which 
varies very slightly with p and, therefore, the recovery factor 
Tg is also seen to be almost constant for the v/hole range of p, 
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Piirther, the total recovery factor will increase in the 
presence of the positive displacement speed, since both the 
first and second-order recovery factors have the same sign* 

The second-order contriVxtion of' cxtomal vgrticit? to 
the heat transfer Gr^(O) is shown in Figure 3.10, It is 
interesting to compare these results with those of Afzal 
and Oberai (1972). For A = 0 and 1, Afzal and Oberai (1972) 
have shown that the heat transfer due to vorticity for a 
given Pr contains an infinite discontinuity for some cri- 
tical value of p. However, in the present work, for the 
case of A= 2 §, we do not observe any such infinite discon- 
tinuity, Thus, in contrast to the results of Afzal and 
Oberai (1972), we note that if we consider- the case of the 
prescribed wall temperature such that full similarity 
(with dissipation included) is preserved (i.e, A= 2p) , then 
the heat transfer G|(0) due to vorticity is well-behaved. 
Figure 5-.. fD shows that for a given Pr, the heat transfer 
G^(0) decreases (in magnitude) as p increases. Since the 
sign of the first and second-order contributions to heat 
transfer is s^e, the total heat transfer increases in the 
presence of external vorticity. The second-order recovery 
factor r^ due to external vorticity is shown in Figure 
and we observe that the recovery factor rg is well-behaved 
in the sense that no singularity of the type mentioned earlier 
is encountered. For a given p, the recovery factor r^ 
increases as Pr increases while for a fixed Pr, the recovery 
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factor rg decreases as p increases, since both the first 
and second-order recovery factors are positive, the overall 
recovery factor increases due to external vorticity. 

For the temperature gradient (in on-coming stream) 
problem, the heat transfer S^(o) is same as reported by 
Afsal and Oberal (1972) while is identically sero. The 
second-order recovery factor r^ due to temperature gradient 
is Shown in Plgure 3.12 . As Pr (alternately, p) increases 
for a given p (given Pr) , the recovery factor r^ decreases. 
The sign of the first and second-order recovery factors being 
same, we conclude that the overall recovery factor will 

increase in the presence of temperature gradient in the 
on-'-comlng stream. 


In the above discussion, we considered the individual 
second-order effects. We shall now show a comparison of the 
second-order results with the first-order in order to indicate 
the relative Importance of the second-order solutions. Per 
the eiample considered here, we choose the same values of the 
various parameters as will be done for the case of gases in 
^3.11, so that the two oases of llguld and gas can also be 
compared with each other. We choose p = 0.5, Pr =0.7, 

Ai-At and A 4 =0 (jointly self-similar case): the corres- 
ponding heat transfer coefficient, Cj^, defined hy (3. 42 ), is 
given by 


Cl. = 


0.6665 - .2716 + 0.5683 B^ R”'^/^ + 0.333 2 

+ .0.5023 B R~V2 ^ oCR""*' ' 
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If we choose etc. to be 0,1, we find that 

the second-order contributions due to longitudinal curva- 
ture is about 4% > transverse curvature 8*^ , displacement 
speed 5% and external vorticity 7*5% when compared to the 
first-order contribution, further, as parameters B^R”^'^^ 

etc. increase , .the relative magnitude of the second-order 

— 1/2 

effects increases; when B^R ' etc. become of the order 
unity, the expression (3.45) shows that the second-order 
contributions are as important as the first-order itself 
and as such the second-order theory dose not hold good. 
Therefore, the above results may be used only when the para- 
meters etc are less than order unity. The corres- 

ponding result for the recovery factor gives the following 
expression: 

r^ = 0.8356-3.840 B^R"**/^ - .0758 B^R"*”^/^ + 1 ,6713 " 

+ 2.594 B^R“^/2 + 1.789 -t-0(R“'’) (3.46) 

Here, as compared to the first-order, the contribution due 
to longitudinal curvature is 45%, transverse curvature 1%, 
displacement speed 4% , external .vorticity 21% and tem- 
perature gradient (in oncoming stream) 21%, It may be seen 
from equation (3.46) that we have to restrict B^R’"^'^^ etc. 
to values less than 0.2; thus for an insulated wall the 
pareimeters etc. are restricted to a much lower 

value as compared to the case of a prescribed wall tempera- 
ture. However, these comparisons are for the particular case 
considered here and no general conclusions may be drawn 
from the limited discussion given above. 
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^^5 II HEAT TRANSFER IN G-ASTC R 

Bopyidary-lBYer Eomtim is for Self-Sim i.l.nr -Solutions 

The thermal bound. .ry-layer equations for gases are 
given by (2.33) and (2,34) as listed in §2.5 of Chapter 2. The 
procedure of obtaining the equations for self-similar solutions 
IS essentially similar to that described for liquids in Part I. 
If/e, therefore, mention here only the important steps leading 
to the final set of ordinary differential equations. The' 

first and second-order temperatures are assumed to be of the 
following form; 


t^+u^/2 - {.TV-H.^( 0 )] g(cr,T7) +H^(0), (3.47) 

^2 ^^1,^2 = H(cr,Ty). (3,48) 

Substituting the above expressions in (2.33) and (2.34), the 

governing thermal boundary- layer equations are reduced to the 
following forms; 


Pir St -Order; 


Pr g + f g I'g- 2(r (f ’ij^-f^’) = ]i;(Pr"'t.l)(f 'f ^ 

g(cr,0)=1, g((r,°o)=o. (3.49) 

Second-Order; 

Er- ■'5 " + f e ' -A f '5 _ 2cr ( f '5^ _ f ’ ) _ A 

• ’ I 

-Pg' +20- (P ’icr-Pcr g') -B(Pr'-''-l)(f'F»)" 
-Bt[Pr’’''(?;i') ' +B(.Pr~1-1)^7](f ’f ") ' + 3 f 'f "}“J 

&((r, 0) s= 0, Bjj('?]-«){) as T/ oo (3.50) 


VH1 
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Here Aj _jj » Bjl ® same as defined in Part I 

w 

and no further comment is needed. 

To obtain the self-similar equations, we assume that 
g and G are of the same form as was assumed for liquids, viz, 
i((r,7?) = 1^(7?) +Ei 2 (T?), (3.51) 

■ 5(cr,7]) = 2 +Gp'^^(7?)] , i = l,t,d,v,H (3.52) 

Substituting the above forms of g and G in (3.49) and (3.50), 
we obtain the following system of ordinary differential equa- 
tions. for the first-order and the second-order heat transfer 
for gases: 

Pirst-Order ; 

Pr"^i^' + fg' -Af =0, 

£-1(0) = 1, i-i(«>) = 0. (3.53) 

Pr~'*i2 + fig “ Spf'gg = (Pr"'’-l)(f ‘f 'V , 

. ggCO) =0 =g 2 (°°). (3.54) 

Second-Order: 

Pr''^‘&”(i)+fG;^i)-(AVvi)f’^^^^+ (1+Ai)p(i)g^'-AP'^i^g^ = Qi, 

G^^\o) = 0, G^^^(7?) = b^ asT]--. (3.55) 

Pr-'^lGg^^U “ (2p+Ai) + (l+Ai)P^i^g2 

_2^'(i)g2_ (Pr--'‘-1)(f 'p’(i))"= Hi, (3.56) 

02(0) =0 =G2(“), (3.56) 

where Xi> Qj,> ^i various second-order effects are 

given in Table 3.3i 
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• Table 5.3s The values ofA’s, Q’s, b*s and R*s in the 
second-order thermal boundary-layer 
equations for gases. 


Second-order i A Q-j 
effect ^ ^ 




N 


i 


Longitudinal 

curvature 


1 -■ST-Hrji\)’ 0 


-Pr~'*(?)g2) '+(Pr“Vl)f 'f " 


Transverse 

ctirvature 


Displacement 

speed 


t 0 -Pr-’l(77i2)'+(Pr“'l-1) 

[-'>7(f’f’’)’-3f'f"] 
d Ad 0 0 0 


External v l-2p 0 0 0 

vorticity 

External enthal- H1-A 0 -rj -a 0 

py gradient 


Again we remark that the full similarity with dissipation is 
obtained with A = 2p. Further, it may be noted that the ther- 
mal boundary-layer equations corresponding to the prescribed 
wall temperature (i.e, for g^ and are the same as 

obtained earlier for liquids (i.e. for g^ and &^^^)in Part I, 
thus the solution for is the same as that for given 

in Part I. 

5.9 Heat Transfer and Recovery Factor ? 

Following the same procedure as in §3.5 for liquids, 
we write down the expression for total heat transfer at the 
wall. Thus, 
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(1 U^(s,0)R-V2|'T^_jj^(0)j(2cr)-V2 

[i^(0) -hSg'^Co) +R-'’/2 2 G:‘(i)(o) +BG2^i)(0)}] , 

i = l,t,d,v,H (3.57) 

Putting q « 0 in (5.57) and using (5.36), we obtain the 
iollowing Gxpnession fon the necoveny f actons j 

=s= 1-2g2(0)/g^(0) (3.58) 

2ti2(0) e;(O(0)-i;(0) &2<i\o)]/i;2(o), 

i = l,t,d,v (3.59a) 

^2 2S.|(0)/gJ(0). (3.59b) 

The heat transfer ooefflolent, with A = 2p, is defined as 
Ci,= (l(2<r)''/2pr Rl/2[ran.|(s,0) (T^Tj,)] 

= -i;(0) - R-V2 2 s;(i)(0) (3.60) 

i=l , . . V ' 

Mact Solution f or the Displacement Sneed Problem ; 

As in the case of liquids, here again we can express 
the solution of the displacement speed problem forA^j = 0 in ' 
terms of the first-order solution. Thus, 
pCd) ^ 

= 7)i;/2, (3.61) 

= 4ig +'!?i^2. 

Using the above solution,, wo pan again express the second-order 

recovery factor in teims of the first~order recovery factor as 
r(d) _ ^ 

2 ^ { 3 » 62 ) 
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3.11 Discussion : 

The heat transfer at the wall and recovery factor are 
obtained in the same way as described in §3.7 of Part I, As 
mentioned earlier, the solution for is identical with 
that of obtained for liquids and, therefore, it is unneces- 
sary to discuss it here. We shall present results for the 
second-order recovery factor r^ only and discuss the contri- 
butions of individual second-order effects. 

The contribution of the longitudinal curvature to the 
recovery factor r^ is shown in Pigure 3.13. I'or a given Pr 
and A]|_ > Of vie observe that the recovery factor ^2 decreases 
(in magnitude) as p increases. Purther, forAq = 0 and 
Pr > .7, the second-order recovery factor T 2 is seen to be 
negative and its sign is, therefore, opposite to that of the 
first-order recovery factor and we conclude that in presence 
of the convex longitudinal curvature, the total recovery 
factor decreases as Pr increases. 

The second-order effect of transverse curvature on 
the recovery factor T 2 is shown in Pigure 3.14. fixed 
At > 0, the recovery factor r^ is positive for Pr > 1 and . 
negative for Pr < 1. Thus we conclude that positive transverse 
curvature will increase the overall recovery factor r^ for 
Pr > 1 and decrease for Pr < 1, 

The effect of displacement on the second-order recovery 
factor rg is shown in Pigure 3.15:. Per A^ 2 recovery 

factor ^2 is positive for Pr > 1 and negative for Pr < 1 , 
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Since the first-order recovery factor is positive, the 
second-order recovery factor rg will add to the first-order 
recovery factor r^ for Pr > 1 whereas the total recovery 
factor r^ will decrease for Pr < 1. 

The second-order contribution due to external vorticity 
to.recovery factor r2, as function of p with Pr as the parameter 
is shown in Figure 3.16. For a given Pr, the second-order 
contributions are larger for negative p as compared to those 
for positive p. The recovery factor is positive for 
Pr > 1 and negative for Pr < 1, As compared to the first-order 
recovery factor r^, the sign of the second-order recovery 
factor r^ is same for Pr > 1 and hence, the total recovery 
factor r^ increases (for Pr < 1, the total recovery factor r^ 
decreases) , 

For the stagnation enthalpy problem, the second-order 
recovery factor r^ is identical to r^, obtained for liquids 
in Part I (shown in Figure 3*12);no further comments are 
needed. 

As in Part I, we present now a comparison of the first- 
order with the second-order solutions. We consider a practical 
case of stagnation region of an axisymmetric body flying in 
air. Thus, p = 0.5, Fr = andAd = 0 (jointly 

self-similar case). Yfe note that the same set of parameters 
were used for the case of liquids also. Further, the solution 
for being the same as that for discussed in Part I, we ‘ 
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shall present here results for the recovery factor ^2 only. 
iFor the case mentioned above, the total recovery factor r^ 
is given by 

= 0.8356-0.1670 - 0.0758 

- .3286 B^jR-V 2 _ 0.8031 ByR“'‘^/2 

- 0.351 BgR-1/2 + o(R-‘’) (3.63) 

If we again choose etc. to be 0.1, we find that as 

compared to the first-order, the second-order contribution 
due to longitudinal curvature is about 2% , transverse 
curvature 1%, displacement speed 4%, external vorticity 
10% and external enthalpy gradient 4%. Equation (3.63) 
suggests that to obtain meaningful results for the second- 
order effects, B-j^R etc. are to be restricted to values 
less than order unity. 

The main results of Part I and Part II are contained in 
references [3] and [ 4 ] respectively. 



CHAPTER 4 


LOOAIiLY-SIMILAR SOIUTIOHS TO THE SECOND-ORDER BOUNDARY- 
LAYER EQUATIONS POR TWO-DKENSIONAL PLOWS 

4*1 Introd-uction : 

In the previous chapter, we studied the self-similar 
solutions of the boundary-layer equations. In various 
engineering problems, however, there are flow situations 
where the similarity conditions are not met and, therefore, 
one has to deal with the nonsimilar flows. Nevertheless, 
the self-similar solutions are of considerable interest since 
they represent the leading term of the locally-similar 
expansion. In literature various methods are available to 
deal with the nonsimilar flows; a brief discussion of these 
methods, with relative merits and demerits, is given in 
Chapter 1 and as mentioned there, we adopt here the (Rjrtler 
series method to construct the locally-similar solutions of 
the second-order boundary-layer flows. In this method, the 
stream function and temperature are represented by a power 
series having streamwise coordinate cr as the expansion 
parameter. Likewise, the pressure gradient parameter and the 
second-order parameters are also expanded in terms of cr . 

The method is applied to uniform flow past two-dimensional 
bodies where the second-order contributors present are lon- 
gitudinal curvature and displacement speed. Although the 
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treatment of the problem, in its essence, is quite general, 
detailed solutions are presented only for a blunted wedge 
geometry. 

We consider, in detail, the boundary -layer solutions 
for a steady uniform flow (free of vorticity) past a semi- 
infinite, symmetric blunted wedge of Included angle Pbcin: such 
that 0 ^ ^ 1 (see Figure 4.1). Parabola ( p£x> = O) and 

flat plate normal to flow direction (^ = l) are the two 
special degenerate cases of the blunted wedge. The boundary- 
layer flow on a blunted wedge starts from a two-dimensional 
stagnation point, then accelerates in the downstream direc- 
tion and finally approaches that past a wedge of Included 
angle Thus, at the stagnation point (o' = O) and far 

downstream (cr-* 0 °), the flow is similar while in the inter- 
mediate region, it has to be regarded as nonsimilar. The 
first-order velocity boundary-layer on a blunted wedge has 
been studied by Chen et.al, (1969) while the thermal boundary- 
layer was considered by Chen (1970), Further, the leading 
term (corresponding to the stagnation point) of the second- 
order contribution for a parabola has been obtained by 
Van i^ke ■(1964a). 

The Gbrtler series method is used to obtain skin 
friction and heat transfer distribution on a blunted wedge. 

Out of the two second-order effects present, the longitudinal 
curvature effect is local in nature and depends only on the 
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first-order boundary-layer quantities. However, the displace- 
ment effect requires the evaluation of the second-order outer 
velocity thus has an upstream influence; it is 

in this spirit that displacement effect is said to be global 
in nature and exhibits the elliptic nature of the original 
lavier-Stokes equations. Methods for obtaining U2(s, O) are 
discussed in Chapter 1 and as mentioned there, we solve the 
perturbed Euler equation with appropriate boundary conditions 
to get the expression of U^Cs, 0) for a blunted wedge. The 
first five terms of the power series for skin-friction and 
heat transfer are evaluated for both the second-order contribu- 
tors. Since the radius of convergence of the series is 
limited, we Eulerize the series to extend its range of vali- 
dity to infinity. The results for large streamwise distances 
are further Improved by accelerating the convergence of the 
series by using either Shanks non-linear transformation or a 
polynomial fit to the partial sums of the series. The results 
of skin-friction and heat transfer distribution for various 
values of wedge angles are displayed graphically. The skin 
friction distribution so obtained for the case of a parabola 
= O) is compared with that obtained via exact numerical 
solutions of the Navier-Stokes equations, by Davis (1972), 
Dennis and Walsh (1971) and Botta et.al. (1972) and the 
discrepancy is found to be less than two percent for Reynolds 
number as low as 10^. 
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Eqxiatlons of t he Bomi(iarv-.La.Yer Flow i 

Consider an incompressible uniform flow past a semi- 
infinite two-dimensional body with prescribed wall temperature. 
The first-order and the second-order velocity and thermal 
(without dissipation terms) boundary-layer equations may be 

written as follows (notations employed are same as used in 
Chapter 2) : 


FjJSt-Order BoundA ry Laver Squat 
Continuityj u^g + = 0 

Uiu^g + = U^(s,0)U-js(s,0) 


Momentum: 


Energy; 


u- 


Hs ” ^3? 


~1 


=* 0 


Boundary oondltions:ui(g,0) =0 =Vi(s,0), hi(s,0) =T 
Matching conditions:u-^(s,N) = TJ'j(s,o) 

h^(s,N) = Hi(o) 

^cpnd-Order Bound ary-laver EaufltinnR ; 


w 


rr 


Continuity: 

Momentum: 


^2s + “ -K(lTv^)lf 

UiU2s + v^ugu + U 2 U-IS + V2U^ii 

= K[lTuiu^g - NUi(s,0)Uis(s,0) 

+Uiu-u^v^]-[EEnjf(s,0) 


+K / u2(s,0)-u| dN] 


s 


Energy; 


+[u^(s,o)U2(s,o)]g 

^2s ^2N ^2 ^ls '^2 


^2NU = h^g + Pr"'* h^jf) 

Boundary conditions: U2(s,0) «0 =.V2(s,0), h 2 (s,o) =0 

Matching conditions: U 2 (s,N)s- - ia)TU^(s,o) +U 2 (s,o), 

^2(3,1) « H^( 0 )^^(s,]!T) 


'(4.1a) 
(4.1b) 
( 4 . 1c) 
(4.1d)‘ 

( 4 . 1e) 


(4.2a) 


00 


(4.2b) 

(4; 2c) 
(4. 2d) 

(4.2e) 



64 


It may Toe noted that the momentum transfer equations 
are same as those given hy Narasimha and Ojha (1967). 

4.3 Kfonsimilar Solutions to the Momentum Transfer Problem : 

We shall use the Gortler type variables, defined as, 
d- = /®U.(s,0)ds , (2<r)“''/^lT^(s, 0) N. (4.3) 

The first and second-order stream functions and 
introduced so that the equations of continuity (4. la) and 


(4.2a) are automatically satisfied. Thus, 

Hi = , T'1=-1/^1S 

H2 = > ■’’g + = -*fes (4.5) 

We assume the following forms for the first and second-order 
stream functions 

= ( 2 T)V 2 f(cr,r)) (4.6) 

^2 = (2o-)''/2i'(cr,-?)) (4.7) 


SO that the momentum equations (4.1b) and (4.2b) are reduced 
to the following forms: 

first-order velocity boundary-layer equations: 

f ff ” (l-f‘2) .> 2cr(f‘f^- l^f") =0, (4.8a) 

f((r, 0) +2crf^(cr,0) =0, (4.8b) 

f’(cr,o) =0, f’(<r,~) =1. (4.8c,d) 

HereAu^ is same as defined earlier by (3.5). 

Second-order velocity boundary-layer equations: 
f’"+ ff" -2Au^ f’f' -f f'f' -2^r (f^’-f + f’f^- f^f”) 

= Bi[-7?(l4Au.j)f*"+(AB^+Au^-1)(f” + ff’+ 2crf’f^ 

+ 4cr ^ f f^dl?) (2 Au.^+ABi)(%.j + a + 20 " a^)]/(1+Auj> 
~ ®d^ * (4. 9a) 
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pCct, 0) + 2cr (cr, 0) = 0, (4.9b) 

P-'Ccr, 0) = 0, F’(<r, ri)^ - B^?)+ asT?- oo (4.9c,d) 

Here the quantities and B^ are given by, respectively 

(5.9b) and (3.9d) with 3 =0 while defined 

by (3.10)« It may be emphasized again that the generalized 

local pressijre gradient parameter An > longitudinal curva- 

ture» parameter At 3 » displacement speed parameter An 

-^d 

are functions of streamwise coordinate cr and the form of 
these parameters varies for different flov/ conditions and 
geometric configurations. We assume the following general 
form of these parameters 




o-P” 

(4.10a) 

^1 

m=0 



Abi 

= 2 A-, 

m=0 


(4.10b) 

A-r^ 

= 2 


(4.10c) 


m=0 '"m 

where p, q and 1 are positive or negative numbers. lo study 
the locally-slmilar solutions of the second-order problem 
(2.8)-X2.9), v/e assume, in addition to (4. 10), the following 
power series for f(cr,'^) and 'Eicr,7)): 


t(<r,ri) = z y>j) 0- 


pm 


m=0 


(4.11) 


F(cr,7|) 


+ ®d crP'^(4.12) 


m=0 


Substituting the above expansions (4, lo)-(4. 1 2) in 
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equations (4.8)-(4#9) and collecting coefficients of the 
like powers of cr, we obtain the equations for the successive 
approximations. It may seem desirable to reduce the coeffi- 
cient functions lAs » P-i » Pp* • • 5 A -i > Ai 

Ad ^ Ad ***^ linear combination of (universal) functions 

that do not depend on (i = 1>2,.,,n), However, 

^ i i 

it was realized that the velocity boundary-layer alone gives 
rise to 89 equations as against the original 16 equations 
while for the thermal bo\indary-layer ; the number of equations 
becomes even Itrger. It woiild be a formidable job to handle 
such a large n-umber of equations and tabulate the results, . 

Here we shall be content with solving the fewer original 
equations for the problem at hand and.no attempt will be 
made to tabulate the universal functions. Instead, we illus- 
trate the procedure in detail for a particular problem, des- 
cribed in the next section. 

4*4 Velocity Bo under v-Laver on a Blunted Wedge ; 

¥e describe here the procedure for studying the locally- 
similar solutions of the velocity boundary-layer due to uniform 
(zero vorticity), steady, laminar flow past a blunted wedge 
shovm in figure 4.1. Het the wedge half-angle be 
(0< poo<1) and its nose radius unity. In the physical plane 
Z = x+iy, the coordinates (s,n) are used to denote respectively 
distances along the surface of the body and normal to it. 
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4.4«1 girst~Order Prolplem : 

First of all, we have to obtain solutions to the first- 
order outer problem; in particular we are interested in 
inviscid flow velocity, U^(s,0).For this purpose, conformal 
mapping is used which transfoims a sharp wedge in the physical 
plane Z = x+iy into the imaginary axis X = 0 of a plane 
C=X+iS- The transformation is 

2 = (1- fe.) P~)- (4.13) 

As shown in Figure 4.1, this transformation maps the blunted 
wedge of include angle into the line 1* length 
element ds = (dx^+dy^)^/^ in Z-plane is mapped to 


dh = in -plane by 

i = Ifll = (4.14) 

SO that along the surface of the blunted v/edge (X= 1) we have 

p. (1--PJ/2 

ds = d4 . (4.15) 

The complex flov/ potential is given by 

<|> = - -1 ( 4 , 16 ) 

from which we obtain the following expressions for the real 
velocity potential, (f) and the stream function, ^1 

(i> = - I (1-pJ[(>i-"l)^- (4.17) 

% = - (1-PJ<(^-1) (4-18) 


The first-order outer velocity along the surface of the blunted 
wedge, as obtained from (4,17), is given by 
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U^(s,0) 


Fs 


%=^ 


W*ds 


1 = C(1+V) 

x=i 


•(1-ftJ/2 


(+.19) 


In "the G-ortler method, since the real potential along the 
surface of the body is used as the streamwise coordinate 
(denoted by cr), we obtain or by evaluating (4.17) at OC = 1; thus, 

= (1-pJ^ ^2. (4.20) 


The generalized local pressure gradient parameter , An 
^ 1 
blunted wedge is now obtained by using (5.5); we have 


Au H [2(r/u.(s,.o)] du.(o-, o)/dcr 

1 ' 

= (l-p^+2p^cr) (l-p^+2(y)“'' (4.21) 

which can be expanded for small cT to yield 

Au ^ ^ pL0-®= p + S (•-2)“(l-pj'‘‘“^V”^ with p. =1. (4.22) 

1 m=0 m=1 . 

The solution of the first-order boundary-layer equation (4.8) 
is expressed as a power series in cT with coefficients depend- 
ing on 17 only; thus, 


f(cr,7?) = 2 f (77) cr“ (4.23) 

m=0 “ 

Substituting (4.22) and (4*23) in (4.8) and collecting coeffi- 
cients of the like pov/ers of o", we obtain the follov(?ing recur- 
sive system of ordinary differential equations; 


Iff ** y ’ 

ftf y V • * y ^ ” 


m 




(4. 24a) 


m-1 m-j 
+ E E 
j=1 i=1 


t t t'm-l 




m-1 
+ E 
D=1 


j_pQ+2(m-j)j 


t I 

1 m— T 


m-1 „ 

- E (l+2j)f,f„ 


(m = 1 ,2. '5.. . . 


^ (a 0/l•K^ 
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with the hoitndary conditions, 

fo(0) = 1^(0) = 0, = 1, 

=0’ =0 (m=1,2,3... .) (4.24o) 

Here e<iuation (4.24a) is the well-lcnown Palkner-Skan equation 
and except for it, all the other equations of the above system 
are -linear. It may he noted here that for = 1 , the expan- 
sion (4.21) reduces toA^ = 1 so that all pj^s (m^l) are zero; 
consequently the expansion of f((r,7)) will have only the lead- 
ing teimi ^qC”^) and the other terms ^^ 01 ) m>1 are absent, 

This flow situation corresponds to the stagnation flow (Hiemenz 
flow) v?hose solution may be obtained from (4.24a). 

The solution for 0 _< < 1 is obtained by solving the 

set of equations (4.24). Gortler (1957) has tabulated the 
results for the first six f’s in (4.23), by splitting them 
into linear combination of universal functions so as to make 
the solution independent of p^s (m21) and depend on p^ only; 
the tabulated results are available for Pq = 0 and 1. ]?or the 
present case, we are interested in the solutions for Pq = 1 
only. With p’s known from (4.22), the first-order solution 
is known in the series form (4.23) upto first six terms. This 
solution is needed to evaluate the displacement effect, 

4.4.2 Second-Order Problem ? 

Of the two second-order effects present, we first consider 
that due to the displacement. In order to study the displacement 
speed problem, we need to know the second-order outer velocity, 
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To evaluate it, we have to first solve the second-order outer 
flow problem (2.28) which may be expressed in the transformed 
coordinates ) as 

(3^/ + (3^/a^^) = 0 {4.25a) 

with the boundary conditions 

1) = 11m § lijn [f(e^,n)-4] (4.25b) 

U oo * Vj->oo 

’vI /2 = ^ upstream. 

Here eq.uation (4.25a) is the Laplace equation in transformed 
coordinates, and the boundary condition (4.25b) specifies that 
the value of '^2 surface of the blimted wedge ('X= 1) is 

governed by the first-order boundary-layer quantities while 
the boundary conditions (4.25c) implies that the second-order 
outer velocity vanishes far upstream. Before proceeding to 
solve the above equation, we evaluate the right hand side of 
the boundary condition (4.25b) as follows: 

Using (4.23), we rewrite the boundary condition (4.25b) as, 

where = lim |')7-f^('J7) | 

Here f^(T?) as ■>7'*°° are obtained from the first-order boundary- 
layer solutions. Using the asymptotic forms of the universal 
functions along with ^’s given by (4.22), we get 

f^(oa) = -.367830, 

= -.276744 + .633740 (l-pj“^ 
f^(oo) =: -.227944 + .929008 d-pj"*’ - 1.127664 (1-pJ“^, 
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f^(o*) = -.204192 + 1.147280 (1-pJ“'' ^2.414976 
+2.047952 

f^(co) = -.196640 + 1.381760 (l-pj"*' -3.956864 (1-pJ"^ 
+5.703264 (1~PJ“^ -3.772288 (l-pj""^, 

= .647901. (4.27) 

Substitution of (4.27) in (4*26) yields tbe expression for 

1) in a series form- whose coefficients are multiplied, by 
§2®+1(ni = 0,1 j 2,3,. . .,). However, it is observed that the series 
converges foi § less than unity due to a singularity at 
+ i. It is possible to enlarge the radius of convergence 
by applying Euler transformation (see §4.5 for detailed discussion) 

p p 

The series (4.26) is recast in powers of a new variable z=% ^(1+^ ) 
Since '^ 2 (C > function of ^ whereas z is even, we 

extract § from the series before introducing the z variable. The 
resulting series for ''p' 2 (^» 1 ) in terms of z variable is given by 

'I) = 1*647901 +.183915 (1-pJ z + [,025480 (1-pJ 

+.069186 (1-pJ^Jz^ + [.008003 (l-pj + .022246 (l-pj^ 
+.028493 (1-Pj^]z^ + [.003487 (l-pj + .010116 (l-pj^ 
+.013774 (1~PJ^ + .01272 d-pj^^z^ + t.001819d-pj 
+.005505 (1-pJ^ + .007790(l-pJ^ + .007868 (l-pj"^ 
+.006145 (1-pJ^z^ + ...} (4.28) 

How we shall solve the second-order outer flow problem 
governed by (4.25a) with boundary conditions (4.25c) and (4.28). 
Instead of solving the equation formally, we shall write down 
the solution for'^ 2 (^»^) inspection. The terms in (4.28) 
are replaced by functions of ^ and X so that they satisfy 
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equation (4*25a) and the boundary conditions (4.25c-) and 
(4.28), In order to satisfy the laplace's equation (4,25a), 
all the terms in '4'g(f,^0 must be harmonic functions and further 
they should reduce to the form (4,28) at X = 1 and also have 
proper behaviour far upstream. The first term in (4.28) is 
which is already harmonic and satisfies the boundary conditions 
on the body surface and upstream; we retain it as such. Other 
terms to (4.28) are, however, not harmonic. These are replaced 
by proper harmonic functions which, at X = 1 should reduce to 

the original terms of equation (4.28). The final solution, thus 
obtained, is given by, 

♦ 2 ^?^= .183915 (1_PJ K-f/(g2+X®)] 

+ [.025480 (1-pJ + .069186 (1-Pj2] [g _ 

+5)(/(g -f [.008003 (1-p^ + .022246 (l-p^J^ 

+ .028493 d-pj^] [| - 3§/(g%X®) + 11?X/4(|2 +x2)2 
+ S^-3SxV4(5^+x 2)5] + [.003487 (1-pJ + .010116 (l-pj2 
+ .013774 (l-pj5 + .01272 (l-pj4][5. 4f/(52+^2) 

+ 41&(/8(g2 4X^)2 + 7(g2_5|;(2)/8(j2+;,,2)3 
+JX(g2.x2)/2(52+x2)'''] + [.001819 (1-pJ + .005505 (1-B )Z 

I QO^ 

+ .007790 (1-PJ2 + .007868 (1-pJ'^ + .006145 (l-pj^] 

[?- 5|/(§2+x2) ♦ 5155X/64(f2+x2)2 + 1 25(5 3.^°°2)/ 

64(<^+X®)^ - 17?l(f2_x2)/8(?2+X^)'’- 

- 5?X'‘)/16(?2+x,2)5]+...J 

Stow we can obtain the displacement speed 02(2,0) from (4.29). 

Thus , 
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Vs.O) = 


= -VCi-Pj 


2, ( P„-5)/2 


where 


1 

[%+a^z *e2^^*ei^z^*a^z‘^+. . .J 


% = .402486 + .095203 (l_pj .054338 (i.pj 


8i 


rOQ-r 

+ .012279 (1-PJ^+ .003360 (1-pJ^, 

^ ~ .125430 + .343291 (1-pJ + .087384 (l-pj^ 
+ .028978 (l-pj^+ .007681 (l-p 
ag = .062571 + .176210 (1-pJ + . gggeo? (1-pj2 

+ .055989 (1-pJ^+ .013826 (l-pj'*, 
^5 = .035172 + .102948 (l-pj * .141353 (i.p^s’ 

+ .133568 (1-pJ^+ .02458 (l-pj'*', 
^4 = .01819 + .005505 (l-pj . .0779 {l-pj2 
+ .07868 (1-pj3 .06145 (l-pj4. 


(4.30) 

The series m the square bracket is divergent for f > 1 (2 > .X) ^ 

To ensure that the series is convergent for the whole range of 

i.e. upto 2 = 1), different factors of the type (l 

(where m was varied) were extracted and convergence wae tested 

for each case; these observations suggested that extracting 
out a factor g(i^2) ( fe.2)/2 ^ 

senes for the complete range of p„ (o < p„ < 1). it 
noted that we have excluded the case p.. = 1, elnoe in this 
case the boundary condition (4.28) reduces to 1jig(g,i) .. o so 
that the laplace's equation is identically satisfied by 
0 leading to U^CsjO) = o. 
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Using ihe expression (4,29) of n /’c! t / 

U2(s, 0) along with (4.19) 

for U,(s,o).we obtata the following expression for 
\ ( =U2(s, 0)/U^(s, 0))i 

= -«i-u (i*e=^)-V2 ^ ^ 

+ ye). . (aj _ 3a^2 + 3a/8 ^ a„/16).3 ^ (a _ 3^3/2 


I ' 

+ 3a2/8 4 a^/16 + 3aQ/l28)z^ + • • •] 


I-inally, the dlsplaoement speed parameter, A „ 

Br 


~ (2cr/B^) (dB^/dcr ) 


°d ' ” ' d' ' “•^d'' 
is evaluated and written as 


(4.31) 
defined by 

(4.32) 


h <5-2 + A4 a -3 + A 4 


A 

With 

Ad^ = 4(l-pJ-‘' 


cr 




Adg = {2S2 - e 2 ) 


Adj = [3E3 - 23^E^ + E^(j,2 _ ^ 

Ad^ = 3^(1-pJ-'^ [4E^ _ + SEg (e 2 - Bj) 


+ (2 E,jE 2 


E 


1 - %)] 


where 

- ^<,V\, 

®2 “ - 3a, + 3a„)/ag 

E3 = (aj - 432 + 6a, - 

■ = <^4 - 5^2 * lOag - 10a, . 5a„)/a„. (^.33) 

Ihe other second-order effect present Is that due to the 
longrtudlnal curvature and will be considered now. Ihe expression 
for the longitudinal curvature K and B, are given by 
K “ cr]-(3-p»)/2 ^ (2<r)-1/2 j-M ( „w 

» 1 vU^(s, 0) (4.34,35) 
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Now the longitudinal curvature parameter A 3^ is given by 
Ab^ = Ai^3-+Ai^ a-2+Ai^ <r5 <r'^ + .... 

= - 4 <r(l-p.„ + 20-)-'' = - 4 (i-pJ-'' s [-Z/d-pJ]”’""' 0® 

m=1 

(4.37) 

Again we note that as expected, for Poo “ ^ expression (4.34) 
gives K = 0, meaning thereby that there is no contribution due 
to the longitudinal curvature for poo = 1* Thus, the above 
results are used for 0 < p^ < 1. 

In order to study the second-order momentum transfer 
problem, we assume, in addition to expansions (4,33) and (4.37), 
the following expansion for F: 

r(o-,’7) = B, E cr “ + S cr” (4.58) 

in=0 ^ m=0 ^ 

In writing dovm, the above expansion, we have decomposed the 
contributions due to longitudinal curvature and displacement 
speed. Substitution of the expansions (4.33), (4.37) and (4.38) 
in equation (4.9) and collection of the coefficients of the like 
powers of cr (after terms multiplied by and have been 
separated) yields the recursive system of governing equations 
for the displacement speed and longitudinal curvature problems. 

The first five equations for both the second-order effects in 
the range 0 £ p^^^ < 1 are written as follows (superscripts 1 and d 
are suppressed) : 

Displacement Speed: 

rf! f? f I ff 

^0 + ^0^0 -2^0 Vo + ^0 ^0 = -2- 
P„(0) = 0 = F^(0), F^(~) = 1. 


(4.39a) 
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- SPo fX + tX * [(2i+1)(FX' 1 f ■?" •) 

i=at1 1 m-i 1 

— ^ OP -n ^ in— i. 


- (2p„+21)fX_i-2ffX.l] + ? lA^ fX 1 • 

i=1 j=0 3 ^-1-3 

- (2p.+A^ ) f'l' . .1 = _ A Pft 

^iij(O) = 0 = Fjjj(o), =0 (m = 1,2,3,4). 

Iiongitudinal Ciirvature; 


(4.39b) 


"o ^ Vc'-2Po <K*Vo 
^o(O) = 0 = ^(0), P;(^)-.-7J ae 77- «. . 

V™ - 2Po ^X ■ f>o * [{2i-M)(PX' . + f-3?" .) 

i=1 1 Di-i 1 

~(2B +2i) f*P* Pl5’'f’ 1 a. V PA " 

° '-A-i """iVJ Xfijlo 

-(2Pi.Ai^) 

a. ? ^ /4>" m-1 

X i (2i a, o„_. 


(4.40a) 


• ' ^ m m-i , 

°i VP + ^ 2 c. f . f . . 

' ^ i=1 3=0 1 0 m~i-3 

m-l m~i 

+ 2 2 23 c, f. f . . 

ixl 3=1 13 m-i-3 


Vo)=o=rm(o), fX)=o 

= (P„ -3)/(1-pJ, 


(m = 1,2, 3, 4) (4.40b) 


“2 = [(fto-3)/2(1-0] &iH/(l-gJ], 

<=3 = [P»-3)/2(1-pJ] [:-Pg+pf/2 + 2p,/(l-R..)+ e/(l-pj2], 
°4 = CP~-3)/2(1-PJJ C-P3+P^Pg - pX + (2p2-p^)/(l-pJ 
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”1 = ? - 2d^], 

^2-12 [-Sdg + (IPg-p^) dj, 

’’S = 8Pl«2 + 2(6P2-P?)d, + 24p2)dJ, 

where , 

«3 = fj + (23*1) aj *1?p. + ^ 

" ->7^1 ("> = 1,2,3,4). 

— d'ty l^^iiTidgr^layer f^li^T 'acteri sties ; 

In order to study aod desoribo the detailed behaviour of 
the boundary-layer, it Is convenient to define few boundary- 
layer oharacteristios. She wantltlee of particular Interest 

loh have to be determined are the velocity profile, displace- 

mert thickness an(J skin fr lot inn m 

in iriction. We now obtain expressions for 

these quantities for both the first ,:j 4 .u 

-cne ii^st-order and the second-order 

boundary-layer flov/s on a blunted wedge. 

Velocity profile; 

The first-order dimensionless velocity profile is given by 


f (d-.i?) 


5: iJV) <r“. 

m=0 “ 


(4*42) 


The second-order velocity profiles corresponding to the dis- 
placement effect and longitudinal curvature effect' are given by 

= S. z o-m , 

'*m=0 ® (4.43a) 


P ^^^(<r,77) ,^5^5; j'(l)(7jj 

m=0 ^ 


(4.43b) 
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A*- / 

^ - ( 2 <r) D^(a, 0 ) A= - 6 * 2 / 2 ) ^ B^(^^. 8 *) ( 4 .^^) 


* 

1 




(4.47a) 

(4.47b) 

(4.48) 


where 

Ai-6 ^ I dv? - 5*^/2 5^ 

^ 00 0 

Aa-8 = / [1-P ('1 )(v 7)3 d*) - 6*?a* 

We assume the following power series "for 

-m ” 

Usirg expressions (4.47) and (4.48) in conjunction with (4.42) 
and (4.43), and collecting the coefficients of the like powers 

of (T , we obtain the following expressions for the coefficient 
functions of the power series of A*^ 

Longitudinal Curvature; 

= 8*2/2 _ P,2/2 + ^ 


= ®o®* - 

^Ip “ ®o*2 * - llm P 

^ // -♦ 00 ^ 


A 


*e * o *» * 


1, = «o'83 * ®?2 - 


^1, = * 5*83 + 8*2/2 _ jtoi pO){ 7 ?). (4.49) 


V- 


Displacement speed J 


K (’?)]. “ = 0,1, 2, 3.^. 


'8m " rf. 
Skin friction: 


(4.50) 


Ihe total skin friction at the wall, T^, is gn,an 
the following expression; 
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Tw = ••• . ('^•51) 

Substituting the inner expansion (2.13) in (4.51) and. using 
(4.4) to ( 4 . 7 ), ( 4 . 11 ) and (4.12), we obtain the following 
expression for 

r = 2(2c*)“''/^uf(s, 0) S f''(0)(T°^ 

' ^ m=0 “ 

+ E-''/2 5 : .(4.52) 


Eor convenience, we define the skin friction coefficient by 

0j= R''/2(2CT-)''/2Ty2u2(s,0) = T., + E-''/2 + ... . (4.53) 

where Tj^ and T 2 are the first and second-order skin friction 
coefficients. Using (4.52) and (4.53), the expressions for T.^ 
and Tg are readily obtained and may be vnritten as 


with 


= E fj0)<r 


m 


in=0 




Tg .= 41 ) + 


- B 

'2 -“1 


2 P. 
m=0 


'KD 


m 


r: 


(d) 


= B , E P 
^ m=0 


”(d) 


m 


(o)(r“^ 

( 0 ) 0 -“ 


(4.54) 


(4. 55a) 
(4.55b) 


Here and represent the overall contribution of the 

longitudinal curvature and displacement speed to the skin 
friction coefficient C^. Por the blunted wedge, using expressions 
( 4 . 31 ) and ( 4 . 35 ), we may express the qtuantities and by 
a power series in <r . Thus, 

Bj = 2 B, <r“ = E.a- 

m=0 m 

+ BjO-^ + 16 d-pJ-''l)^<r ++...] 

(4.56) 
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^ mfo \i <r I 

+ 4(1-0"^ 0- ^-8(l-ft„)-V5 + 16(1-Rj-4ff-4^.__| 

Introducing the above expressiona in (4.45), we obtain the 
following expressions for andT^^^* 

4* 2 ^ 


T 


( 1 ) 


[.*(1) 


in=0 


m 


with P^l) = I B 

n=0 n ® " 

= 2 P*(d) m 

^ m=0 

with = I B 


(4.58) 


(4.59) 


Here B^^'a and B^'a are obtainable from (4.46) and (4.47). 


4.5 Heat tra nsfer Prohl^m . 


In the preceding sections, we considered the nonslmllar 
solutions of the momentum transfer problem (§4.3) and studied 
explicitly the example of a blunted wedge (§4.4). ffe now con- 
sider the corresponding thermal problem here and In the next 
section (§4. 6) we shai study the example of the thermal 
"boundary-layer on a blunted wedge. 

Ihe thermal -boundary-layer equations (4.I0) and (4.2o) 
may be rewritten In (<r,V) variables with the help of the 

following transformations for the first and second-order 
enthalpy: 


*“1 “ «(<’■,’?) +H^(0), 

"a = tV^iCO)] G(er,7;). 
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Using the above expressions along with (4.11) and (4.12), the 
thermal boundary-layer equations are reduced to the following 
form: 

Eirst-Order : 

Pr‘“''g'* + fg‘ - 20- (f‘g^-f^g’) = 0, 

g(<r, 0) = 1, g((r, ~) = 0. (4.62) 

Second-Order ; 

Pr"’''G'’ + fo' - 20- (f'Gjj.- f(^G') = ~l'g'+ 2 (r (p’g^-E^g') 

-Pr-''B^(7?g')', 

G((r,0) =0, G(cr,oo) = o. (4.65) 

In order to study the locally- similar solutions of (4.52) and 
(4.53), we assume the following power series for gicTfTj) and 
G(cr,i7): 

g(<r,7?) = s g (7)) S-I® (4-. 64) 

m=0 “ 

G(<r,l3) = S oO)(T)) + (4.65) 

111=0 

Substituting the above series of g(cr,'^) and G(o“,l7) in equations 
(4.62) and (4.65) and collecting the coefficients of the like 
powers of (y~, we obtain the equations for the successive approxi- 
mations. However, as commented earlier for the momentum transfer 
problem, we shall not attempt to obtain the solutions to these 
equations in terms of the universal functions; instead, the 
procedure is fully illustrated in the next section by considering 
the problem of heat transfer on a blunted wedge* 
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iti.er;iial Boimdarv-LavRT- on a Blimtoil 

The momentum transfer on a blunted wedge was studied in 
44.4; the corresponding heat transfer problem Is considered 
here. Although we shall study the constant wall temperature 
case, the method may be readily extended to the case of a 
vaxiable wall temperature. 

Eirst-order rrnhi e>m ♦ 


Following the procedure outlined in the preceding section, 
we assume the following power series form foTg(cr,y}): 


si<r,r)) = 2 g^(7?)(r“ 

m=0 


(4.66) 


Substituting (4.66) In (4.62) and collecting the coefficients 
of the like powers of o' .we obtain the following recursive 
system of ordinary differential equations: 

•n -1 ” » 

8o + foSo = 0> 


go(0) = 1> g(“) = 0. 

•D —1 " ‘ « ni 

®m * ^0 Sm - 2 “ foSm = -.f/2l+1)f^^_^ 


^(0) - 0, g„(~) = 0 (m = 1 , 2 , 3 ,... J (4.67b) 

Equation (4.67a) will yield the well-known similar solution for 
the two-dimensional problem and has been studied quite extensively 
in literature. Further, for Prandtl number unity, the problem 
(4.67) reduces to the one corresponding to the cross-flow 
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problem studied by Gdrtler (1957) and as such the solution may 
be written directly with the help of the universal functions 
tabulated by Gbrtler (1957), This fact has been used by 
Chen (1970) while studying the first-order heat-transfer on 
a blunted wedge. However, Prandtl number unity is not of much 
interest. For air, Prandtl number is close to 0.7 and we shall 
choose this value for the present study. As such, the equations 
(4.67) has been solved by numerical integration to obtain solu- 
tions needed for the study of the second-order problem. 

4,6.2 Second order 

Por the second-order problem, we assume the following 
expansion for G(cr,’?): 

-r”. (4.68) 

Here we have separated the oontributlons due to displacement 
and longitudinal curvature. Substituting the expression (4.68) 

In eciuatlon (4,65) and collecting coefficients of various powers 
of o- , we obtain the following recursive systems of ordinary 

differential equations for the displacement and longitudinal 
curvature problems: 


Displacement speed: 

f t m-1 , 

^^^o - ^^O * 2i[8A-l * 

- t(2it1)(fiVl -^A-i) (So C± - <Vi)] 


~ 2mf *G 

m 0 m 0 m 


m 


m-1 m-i , 

^ ^ 2 A , (f.G , . F e > 

i=1 d=1 ^i 3 m-i-j 


= 0 = 


(m = 0,1, 2, 3, 4) 
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I<ongitudinal curvature: 

,-1 n” 


4 ^0 


( 4 . 70 ) 


‘'m + - a» < + a, _ 

m-l ^ 

i!, r°l Vl + Si Vj - “ t(2i+1) 

i=1 

®m-i - Ci) +Al (*‘o 4-1 - 4 Vi)] 

m-l m~i 

^ 2 At (f. G V \ 

i~1 3=1 0 m-i-j ” ^3 

e^(o) = 0 = Oj,(“). (m = 0,1, 2, 3, 4). 

We now obtain expressions for two Important thermal 

boundary-layer oharaotorlstlos, namely temperature profile and 
heat transfer at the wall. 

Temperature profile: 

first-order and the second-order temperature profiles 

are given by (4.66) and (4.68) and are repeated here for ready 
reference. 

First-order: 


Bio-, 7)) = S ^(7?)cr"> 

m=0 ^ 

(4.56) 

Second-order: 

. 

G^^^(<r,77) ^22 (j(l)('n) (j-m 

^ m=C| ^ ^ 

(4. 58a) 

G^^)(<r,7;) = B 2: G^^>(7)) cr”' 

in?=0 

(4. 58b) 


Heat transfer: 


The local heat transfer at the wall, is defined 
ci„ = -k (3T/0n) 


as. 


n=0 


(4.71) 
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Substituting the inner expansion (2.15) for T in (4.71) and 
using (4.66) and (4.68) we obtain the following expression 


for 


0) 

L 2 4^o)(r “ + 

111=0 m=0 


(4.72) 


For convenience, we define the heat transfer coefficient by 

Pr (2cr)''/2^^ [a)^-H^(0)]-VUi(s,0) 

^ r-1/2 


= 


q2+ 


(4.73) 


where q^ and qg are, respectively, the first-order and the 
second-arder heat transfer coefficients. Oomparison of (4.72) 
and (4.73) leads to the following expressions for and (I2* 


.m 


2 gjj,(0) cr- , 

in=0 ® 


with 


(1)_ 

2 - 

B. Z 
m=0 

g'(l) 

m 

(<i)_ 

B, Z 

s'(1) 

2 ~ 

^ m=C 

m 


(0) <T 


m 


(0) cr®^. 


(4.74) 

(4. 75a) 
(4.75b) 


Here q^^^ and q^^^ represent the overall contribution of the 
longitudinal curvature and displacement speed to the heat 
transfer coefficient 0^^. Introducing power series for 
(4.56) and B^ (4.57) in (4.75)) we rewrite the expressions for 
q^^^ and q^^^ in the following form: 

m . / V 

. (4.76) 




With S^'> 

= E Bw 

m=0 

m 

m 

1« m-n ' 
n=0 n 

(<i)_ V 

n*(a-)^m 


5 « . 
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4*7 Extending Range of Validity of the Series Solutions in 
the Itownstream Direction ; 

The reinge of validity of the series solutions, studied 
in the preceding sections, is restricted in the downstream 
direction because of the limited radius of convergence of the 
series. We now discuss various methods which have been used 
to extend the range of validity of the solutions in the down- 
stream direction. It is well-laiown that in many applied 
mechanics problems, a power-series expansion has a small radius 
of convergence because of a singularity in the complex plane, 
usually on the negative axis, although the expansion has a 
physical significance only for the real positive values of the 
expansion parameter. In such cases, if we know the location 
of the singularity, it is usually possible to enlarge the 
radius of convergence by transferring the singularity farther 
away, possibly to infinity. Depending on the knowledge of the 
location and nature of the singularity, various transformations 
are available to extract maximum information from a slowly 
convergent or even a divergent series. Bellman (1955) and 
Van IDjyke (I964h) have given several interesting examples of 
such transformations, 

Euler transformation is one of the most powerful and 
commonly employed techniques for enlarging the radius of 
convergence of a power-series expansion. According to this 
method, if there is only one singularity and its location is 
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known, then it is possible to transfer the singularity to infinity 
so that the transformed series has an infinite radius of conver- 
gence. For a power-series having one of the coordinates (say cr) 
as the expansion parameters, it is usually found that its radius 
of convergence is limited because of a singularity on the nega- 
tive axis (say at cT = -p) although from physical point of 
view, only the positive values of (T are relevant. Such artifi- 
cial limitation can be overcome by a simple conformal mapping 
whereby the singularity is shifted to infinity - the techniq.ue 
is known as the Euler transformation which consists of introduc- 
ing a new variable z defined as, 

2 = (r/((r+p) (4.78) 

The transfoi?med series in z has the radius of convergence 
z = 1 which corresponds to cr- oo. Thus, the transformed series 
assures convergence for all values of <r , and the utility of 
series solutions is considerably improved. As a first step we 
have, therefore, to estimate the value of p. 

The power-series expansion in the case of a blunted wedge 
have cr as the expansion parameter. Using (4.20), the expression 
(4.19) for the inviscid first-order velocity U^(s,0)may be 
written in the following forms 

U,(s,0) = ^ (4.79) 

The above expression shows that U.|(s, o) is singular at 
(T = -(1-p^)/2. This suggests that the first-order series- 
solutions will also exhibit a singularity at cr = -(l-p^/2. 
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This expectation is further supported from the trend of the 
ratio of the successive coefficients f^,f2...f^ of the first 
order series solutions (e.g,4.54)> which seems to approach 
constant value equal to for a fixed (cf. §4.8). 

Thus, we obtain the following form of the Euler transformationi 

z = cr/((r+(l-pJ/2) (4.80) 

We suggest that the above form of the Euler transformation is 
suitable for the second-order series solutions for the longitu- 
dinal curvature (4.34) and the displacement speed (4.30); both 
the expressions are singular at cr = (l-p^)/2. Again, the 
ratio of the successive coefficients of the second-order series 
solutions [e.g. , m = 0-4 in (4.58) and (4.59)] is found 

to approach a constant value equal to (l-p^)/2 and, therefore, 
suggest that the radius of convergence is limited by a singularity 
at cr = -(l-p^)/2(see §4.8). have used the same Euler transfor- 
mation (4.80) for both the momentum transfer and heat transfer 
problems since the above comments regarding the radius of con- 
vergence of the series solutions hold good for the heat-transfer 
problem as well. 

The Bulerized series though convergent upto z = 1 (i.e. cT-* ~) , 
the convergence is found to be rather slow leading to considerable 
error for large values of cr . To extract further information from 
such a slowly convergent series, it is sometimes possible to 
accelerate the convergence by applying Shanks non-linear transfer- - 
mation. Shanks (1955) has given several examples to show how a 



89 


lot of additional information can be obtained from a few terms 
of a given series. The Shanks method consists of applying the 
following non-linear transformation to the partial sums of the 
series (cf, e»g. Tan Dyke 1964b): denoting n~th, partial sum 
we write, 

“ ^^n+1 ^n-1 ^ ^n~1 ~^^n^ * (4.81) 

The process may be repeated by using c-j(Sj^)'s as the S^^'s for 
the next iteration. 

Another technique is to fit a polynomial in 1/n (n denotes 
the number of terms) to the last few partial sums of the series 
and thus extrapolate the sum of the series for n - °o. Using 
the last three partial sums of a series containing n number of 
terms, the polynomial fit will predict the sum of the series 
given by the following expression: 

S = (n-2)V2 S ^_2 - (n-1)2 + n^/a S^. (4.82) 

In particular, for n = 6 (for the first-order problem where the 
first six terms of the series are known) 

S = 8S^ - 25 + 18Sg (4.82a) 

and for n - 5 (second-order problem) 

S = 4.5 - 16 +12.5 S^. (4.82b) 

T'or the present case, it will be shown in §4.8 that the latter 
method gives better results as compared to the Shanks non-linear 


transformat ion . 
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^•'7*1 Euler ized .skin- friction and heat transfer ; 

Since the series solutions for the skin friction and heat 
transfer are of the same kind, we illustrate the Eulerization 
procedure for a general case. Consider a pov/er—series expansion 
for some variable 0 to be of the following form: 

0 =00 + 0-1^ + 02^ ^ 0’f’^ 04^^ 05^^ + • • • . ( 4. 83) 

Introducing the variable z (4.80) in the above series and 
denoting the resulting Eulerized series by 0, we obtain the 
following expression for 0 in the operator fom: 

0 = 1* (0ju»z) = 00 + P0^z + (P0^+P^02)z^ + (p0^+2p^02 + p^0j)z^’ 

+(p0^+3p^02+3p^05+P^0^)z^ + (p0^+4p^02+6p^0^+4pV4 

+5p^05) + .... (4.84) 

where p = (l-.pJ/2. 

Using the above result along with the definition of the operator 
1 given by (4.84), we may readily obtain Eulerized series solu- 
tions for the skin-friction and heat transfer. 

Eirst-order r 

Eulerization of (4.54) and (4.74) yields the following 
Eulerized series solutions for the first-order skin friction and 
heat transfer: 

Ti = 1 [f^ (0) , a] 

= iX (0) , a] 


(m = 0-5) 


(4.85) 

( 4 . 86 ) 
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liongitudlnal Gurvattire ? 

Bulerizing (4.58) and (4.76), we obtain the following 
Eulerized series solutions for the overall contribution of the 
longitudinal curvature to slcin friction and heat transfer! 

z] (4.87) 

(m = 0-4) (4.88) 

Disnlacement Speed ? 


Bulerizing (4.59) and (4,77), we obtain the following 
Bulerized series solutions for the overall contribution of the 
displacement speed to skin friction and heat transfer; 

•*'^*(d) 


f(d) 

2 


z] 


(m = 0-4) 


(4.89) 

(4.90) 


To extract further information from the Eulerized series 
solutions, both the methods described earlier, i*e. that of 
Shanks non-linear transformation and extrapolation, were used. 

The results thus improved are described in §4.8 and it will be 
pointed out hov7 much improvement is achieved by the two techni- 
ques of accelerating the convergence of the Eulerized series. 

The total Eulerized skin friction and heat transfer coefficients, 
combined with extrapolation, may be expressed in the following 


Cf = (f(d) +f(l))^ (4.91) 

= q^ + (q^^^ (4.92) 


form: 
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I-Wnerica l Results and Dismissinn i 

The recursive system of differential equations for the 
first and second-order momentum and heat transfer problems are 
solved by numerical integration, using Runge-Kutta-Gill method. 

The solutions are obtained for various values of wedge angles 
( in the range 0 < < 1 . 

We first describe the results for the first-order problem. 
As mentioned earlier (§4.6) Chen (1970) has studied the first- 
order heat transfer for p^ = 0.5 and Pr = i , by utilizing the 
universal functions of the cross flow problem, tabulated by 
Gortler (1957). However, when Pr is different from unity, these 
solutions are to be obtained numerically. In the present calcu- 
lations we take the value of Pr =0.7 as it is the value for air 
and hence of more practical interest, 

Por the first-order boundary-layer, the first six terms in 
momentum and heat transfer are computed from the solutions of 
equations (4.24) and (4.67). It may be noted that the first two 
equations (for both momentum and heat transfer problems) in the 
above system are independent of p^ and hence the solutions of 
these two equations will only be valid in the immediate neigh- 
bourhood o^ the stagnation point. The velocity profiles for three 
representative values of p^, namely p^ =0, 0,5 and 0,8 are 
shown in Pigures 4.2a,b and c. The first six coefficients of 
the skin friction (i.e. fj^'(0), m = 0-5) and displacement 
thickness 6 (i.e. 6 ^, m = 0-5)versus p^ are shown in Pigure 

4.3. The first two coefficients of skin friction fflVo) and flVo)! 
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and displacement thickness (6^ and dg) are independent of (so 
are their governing equations) while others increase monotoni- 
cally as increases. For a fixed value of Fig-ure 4.3 
shows that the ratio of the successive coefficients j 

^n^^n+l ’ ” ~ tends to a constant value, equal to (1-p^)/2 

and supports what was anticipated in 64»7. Thus the series 
solutions (4.45) and (4.54) is expected to have a radius a con- 
vergence of (l-p^/2 and are useful only for cr < (l-p^)/2. 

However, the region of validity of the results is extended in 
the downstream direction by Eulerizing the series with p = (l-p^JI/^. 
If there is no other singularity, in the complex-plane, then the 
Eulerization will lead to results which are valid every where on 
the blunted wedge. To test the Euler ized results for various 
values of p^^, we consider the limit as z 1 (corresponding 
to cr-» oo) where the properties of the boundary- layer should 
approach those of a wedge (given by well-known Ealkner-Skan 
equation with corresponding value of p^^ . The Eulerized series 
of the skin friction (4.85) for z = 1, shows that the successive 
partial s-ums do approach a value equal to that given by Ealkner- 
Skan equation for all p^^ (see Table 4.1). The difference between 
the last partial sum and the exact value is about 15 percent 
for p^ = 0, 0.6 percent for p„ = 0.5 and almost zero for p^ = 1.0. 
To ftirther improve the convergence of the Eulerized series, two 
methods were tried. EirsJ: the non-linear transformation of 
Shanks was applied but the results, although improved, were 
still in error by about 7 percent for poo = 0* ^ polynomial 
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fit in 1/n to the last three partial sums was used and the 
results (cited as improved Eulerized) have shown a marked 
improvement (error is about 1.3 percent for = 0) . A compara- 
tive study of the Bulerized and improved Eulerized results for 
z = 1 with those of Ealkner-Skan solutions of the corresponding 
Pooy obtained by Smith (1954) is shown in Table 4.1 below: 

Table 4*1J The Bulerized first-order skin friction 
results at z = 1. 


Poo 

Eulerized f"^ 

Improved a 
E ulerized 

A.M.O. Smith 
(1954) 

0.0 

0.5410 

0.4636 

0,469600 

0.1 

0.6298 

0.5815 

0.587035 

0.2 

0.7129 

0.6840 

0.686708 

0.3 

0.7909 

0.7743 

0.774755 

0.4 

0.8643 

0.8553 

0.854421 

0.5 

0.9336 

0.9289 

0.927680 

0. 6 

0.9992 

0.9970 

0.995836 

0.7 

1.0616 

1.0606 

1.0598078 

0.8 

1.1211 

1.1207 

1.120267 


The improved Eulerized results for the skin friction versus 
z for various values of are shown in Figure 4.4 and we note 
that the skin friction decreases monotonically for increasing 
streamwise distance for 0 < < 1, and remains constant for 

Poo = 'I- 
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Por the first-order heat transfer problem the temperat-ure 
profile for some representative values of (= 0*5 and 0.8) 

are shown in ^Figures 4.5a,b and c. The results for the coeffi- 
cients of the heat transfer at wall gjjjCo) versus for n = 0-5 

are plotted in Figure 4.6. Except for the first two coefficients, 

! I 

i*e. Sq\0) and g.|(0) which remain constant, other coefficients 
increase monotonically as increases. Further, for a given 
Poo> '^^6 ratio of the successive coefficients, i.e. SjQ(0)/gjjj^^(0) 
for m = 0-5, is found to approach (l-Poo)/2, implying thereby 
that the radius of convergence of the series for heat transfer 
rate is same as that for the skin friction. The Euler ized and 
improved Bulerized results for heat transfer at z = 1 are tested 
and a comparative study is shown in Table 4.2. 

Table 4. 2: Eulerized first-order heat transfer results 


at z = 1. 


Poo 

Eulerized 

Improved 

Eulerized 

Evans ( 1967) 

0.0 

0.4143 

0.4029 

0.41391 

0.1 

0.4275 

0.4238 

0.43140 

0.2 

0.4391 

0.4399 

0.44438 

0.3 

0.4494 

0.4523 

0.45470 

0.4 

0.4585 

0.4621 

0.46323 

0.5 

0.4665 

0.4700 

0.47049 

0.6 

0.4737 

0.4766 

0.47680 

0.7 

0.4801 

0.4823 

0.48235 

0.8 

0.4859 

0.4872 

0.48732 
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Ihe above table shows that for any falue of the Improved 
Eulensed results compare well (within 2 percent) with the 
corresponding solutions for a wedge, given by Evans (l967). 
Figure 4.7 shows the improved Eulerized heat transfer versus 
z for various values of The figure shows that for = 1 , 
the heat transfer is same for all z while for other value of' 

Poo It decreases monotonically as z increases. Further, this 
rate of decrease of heat transfer increases as decreases. 

We now describe the second-order results for the displace- 
ment effect, obtained from the solutions of the first five 
equations in each of the systems (4.59) and (4.69); these results 
are displayed in Figures 4.8 to 4. 14. In Figures 4.8a,b and c 
the coefficients of the power series (4.45a) for the velocity 
profile are presented for = 0, 0.5 and 0.8 respectively. In 
Figure 4.9 the coefficients Fj^’^'i)(o) of the skin friction 
series (4.55) are plotted versus p<„ and the corresponding 
results (A^^ versus pj for the displacement thickness (4,50) 
are shown in Figure 4. 10. lo get a better insight in the struc- 
ture of the displacement speed solutions, the coefficients F*^^^ 
of the power-series (4.59) for the overall displacement speel 
contribution to the skin friction are also plotted as a 

function of p^ in Figure 4.9. It may be seen from Figure 4.9 
that the first coefficient of (4.59) decreases as p^ 

increases; the skin friction at the stagnation point (cr = 0) 
decreases as p^ increases and is zero for p<„ = 1. Further, the 
other coefficients F^j^^^^show a monotonic increase with p^^ 
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Us4ng the series (4.59) and the coefficients from 

m 

Figure 4.9j we can compute the total contribution of the dis- 
placement effect to the skin friction. Further, for a fixed 
Poo> Figure 4.9 shows that the successive ratio 
(m - 0-4) approaches the value (l-ppg)/2» Therefore, for improv- 
ing the convergence of the result s,*Eulerizat ion with p = (l-pJ/2 
is employed. The Eulerized skin friction contribution due to 
the displacement effect f is plotted versus z for various 
values of as is shown in Figure 4.11. For a given p^, the 
second-order contribution to the skin friction is negative near 
the stagnation point (o'‘= z = O) and changes sign for z >0.65 
and becomes zero for z = 1. This behaviour of the skin friction 
is similar to that of the displacement speed. For the particular 
case of Poo = 1 (plane stagnation point flow) the displacement 
speed is zero and so is the contribution to the skin friction etc. 

To study the contribution of displacement to heat transfer, 
the coefficients for the temperature profiles (4.58b) are shown 
in Figures 4.12a,b and c for poo = 0, 0.5 and 0.8. The first 
five coefficients of the power series (4.75), i.e. G-j^^‘^^(0) for 
m = 0-4 versus p^ are shown in Figure 4.13. The overall contri- 
bution of displacement effect to heat transfer is given by 
power series (^.77) whose first five coefficients are also 
displayed in the same figure; the first coefficient of 

series (4.77) decreases as p^ increases, which means that the 
contribution of the displacement speed near cr = 0 (stagnation 
point) decreases as p^ increases, while the other coefficients 
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show an increase as increases. The radius of convergence 
of series (4»77) remains the same as in the momentum transfer 
problem so that the results can be used for cT < (l-Poo)/2* 

Bulerized results for the series (4«77)j given by (4.90), 
are shown in figure 4.14 for various values of The quanti- 
tative behaviour of the heat transfer ^ contribution is 
similar to that of the skin friction and, hence no further 

comments are needed. 

For the longitudinal curvature problem, governed by 
equations (4.40) and (4.70), the results of numerical integration 
of the first five equations for the momentum and heat transfer 
problems are shown in Figures 4.15 -■ 4.21. The velocity profiles 
for = 0, 0.5 and 0.8 are shown in Figures 4.15a, b and c. 

The first five coefficients of the series for the skin friction, 
(i.e. F*'^^^(0) given by(4»55a), and for the displacement thick- 
ness, i.e. given by (4.49), are plotted versus p^o as shown 

in Figures 4.16 and. 4.17 respectively. The overall contribution 
of the longitudinal curvature to the skin friction 7* , given 
by (4.58) is also shown in Figure 4.16 - the first five coeffi- 
cients f”^^^ of the series are plotted as a function of p^^. For 
a given p^, the Figure 4.16 shows thatthe ratio of the successive 

coefficients 

the radius of convergence of this series is the same as that 
mentioned earlier. The Bulerized skin friction due to the longi- 
tudinal curvature , given by (4.87), 1® stio^n in Figure 4.18 

for various values of p^o. a given p<„, this figure shows 
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■fctiat the longitudinal curvature contribution to skin friction 
is negative and decreases (in magnitude) as streamwise distance 
z increases upto O.65. Por z >0.65 the sign becomes positive 
but the magnitude is very small and is zero for z = 1. 

The temperature profile for the second-order contribution 
of longitudinal curvature are plotted in Figures 4.19a,b and c 
for the heat transfer rate at the wall, of the series (4. 75a) 
and (4.76) are shown in Figures 4.20a and 4.20b. The Eulerized 
results for the heat transfer versus z for various values 

of are shown in Figure 4.21. The qualitative behaviour of 
the heat transfer results is the same as that of the skin fric- 
tion and hence no fin?ther comments are needed. 

So far we have discussed the results for each of the 
individual effects. We now describe the Eulerized results for 

A A 

the total skin friction and heat transfer 0^^ and discuss 
their validity with respect to the streamwise coordinate z or cT 
and Reynolds number. For clarity of presentation, the skin 
friction on the blunted wedge for various values of are 
spread over two figures, namely 4.22a and 4.22b and values of 
Reynolds number are 10^, 10^ and 10^, Figure 4.22 shows 
that for a given. p^ the skin friction decreases monotonically 
as z increases for R > 10^. However for lower Reynolds number 
e.g. R = 100 and p > 0.2, the skin friction is found to increase 
initially with increasing z and then decreases to approach the 
Falkner-Skan value at z - 1. A comparison of the present results 
at z = 1, with those of Snith’s Falkner-Skan solutions is given 
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The present skin friction results for a parabola ( Poo = O) 
are compared with the exact numerical solutions of the Navier- 
Stokes equations, obtained independently by Botta et.al.(l97l), 
Ba'vis (1972) and Dennis et.al. (1972). A comparison of the 
present results for the skin friction at the nose of the para- 
bola with those obtained by the above mentioned authors and 
that by Van Dyke (1964) is given in Table 4.3 below: 


Table 4.3: 

Skin friction at the nose of the 

= 0). 

parabola 

= ^/V 

Present 

results 

Davis 

(1972) 

Dennis et.al. 
(1971) 

Botta 
et . al » 

(1972) 

Van Dyke 
(1964) 

10^ 

1.316 

1.427 

1.424 

1.428 

1.314 

10^ 

1.616 

1.626 

1.625 

1.624 

1.607 

10^ 

1.684 

1.710 


1.706 

1.700 

1o5 

1.712 



1.734 

1.730 


Figure 4. 23 shows a comparison between the present results 
and those of Davis (1972) for the distribution of the skin- 
friction over a parabola at various Rejmolds number. Upto 

X 

Reynolds number as low as 10"^, the difference in the two 
results is only about one percent. Even at R = 10^, our 
results are lower by about 7 percent, most of the difference 
in the two being near the stagnation region. Further, Davis 
has also given, v/hat he calls the 'parabolic approximations' 
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and it is pointed out ty Ifevis that the difference between the 
exact results and parabolic approximations is mainly due to the 
displacement effects. We have also shown oirr corresponding 
parabolic approximation results by excluding the contribution 
of the displacement effect to the skin friction. It may be seen 
from Pigure 4*23 that the two results for the 'parabolic approxi- 
mations' compare well (within tv/o percent) even for Reynolds 
number as low as 100. This suggests that we have probably over 
estimated the displacement contribution to the skin friction. 

In Rigure 4.24, we have compared the skin friction results with 
those of Dennis and Welsh (1971) for R = ~ and 100. Also we have 
plotted the results of Smith and Clutter (1963) for R = «>. It 
may be seen that for R = oo, our results are close to those of 
Smith and Clutter whereas the values obtained by Dennis and Ifelsh 
are somev/hat higher for large values of the streamwise distance 
z. Ror R = 100 also our values are lower than those of Dennis 
and Welsh, This behaviour is consistent with that observed 
from Rigure 4.23 (comparison v/ith Davis) as might be expected 
in view of the fact that the results of Davis (1972), Dennis 
et.al. (1972) and Botta et.al. (1971) compare quite well among 
themselves. 

The Eulerized total heat transfer coefficient C for 

Q. 

various values of and R = 10^ and 10^ is plotted 

as a function of the streamwise coordinate z in Rigures 4,25a 
and 4.25b. The effect of variation of and R on the heat 
transfer distribution is qualitatively similar to that described 
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earlier for the corresponding skin friction results and as 
such no further conanent is needed. A test of the present 
results at z = 1 (tr- «>) with those of Evans (1967) for the 
corresponding (see Table 4.2) shows an agreement within 
two percent. 

In view of the above discussion, v\re suggest that the 
present results for the skin friction and heat transfer 
distribution on a blunted wedge may be used v/ith little error 
(within two percent) for Reynolds number upto 10^; for still 
lower Reynolds number, however, one has to be more cautious. 



GHAPIER 5 


NON SHUAR SBOOND-ORNER BOUNDARY-LIYBR SOLUTIONS 
VALID FAR DOWNSTREAM 


5.1 Introduction : 

In the present chapter, we study the locally-slmilar 
solutions valid for large values of streamwise coordinate o' * 
Apart from the fact that these solutions are of interest in 
themselves, they are also useful in assessing the validity of 
the Ghrtler and Blassius type series solutions (obtained for 
small cr ) for large cr * 

Various methods, with varying degree of complexity and 
accuracy, are available for studying the present problem. The 
finite difference and difference differential schemes are 
available to calculate the development of the boundary-layer 
downstream of a specified initial velocity profile. Although 
exact in principal, these methods are unattractive for a 
parametric study due to the large amount of calculations 
involved. An analytical approach is preferred. To obtain 
solutions for cT-* «>, variables are expanded in a complete set 
of eigenfunctions leaving unspecified multiplicative constants, 
which are to be detemined from the prescribed initial condi- 
tions, In principle, this analytical method is exact but in 
practice, the eigenfunction representation is truncated after 
a few terms to yield an approximate solution. This approach 
has been used for the first-order boundary-layer flows by 
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Libby and his several colleagues and a good review has recently 
been made by Libby (1970). In their study, the required set 
of eigenfunctions are determined from the linearized equation 
(obtained by linearizing about the limiting similar solution). 
This linearized equation is, in some sense, similar to the 
well-known Sturm-Liouville equation for which the eigenvalues 
are known to be real and discrete. 


Such a procedure is used here to study the second-order 
momentum transfer problem for large values of o' . We first 
consider the initial value problem for each of the second-order 
effects, namely the longitudinal curvature, transverse curva- 
ture, displacement and external vorticity. We let the pressure- 

gradient parameter An » longitudinal curvature parameter An , 

^1 -"I 

transverse curvature parameter and displacement speed 

t 

parameter An "be constant. Since we are primarily interested 
d 

here in studying the eigenvalue problem, the restriction of 
choosing constant parameters simplifies the problem conside- 
rably without affecting its general nature. At a streamwise 
station an initial velocity profile (upto second-order) 
differing slightly from the limiting similar profile for 
the specified parameters is assumed. We are interested in 
studying the downstream behaviour of the solution so as to 
find how the limiting similar solution is approached as o'"-* «>, 
The results obtained for a variety of parameters are displayed 
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graphically, A critical discussion of the results is given; 
in particular, the effect of parameter variation on the stabi- 
lity of the flow is examined and the implications of the sin- 
gular behaviour of the solutions to the design of body shapes 
having desired aerodynamic characteristics is pointed out. 


The application of the above foimiulation is illustrated 
for a blunted wedge, v/here only the longitudinal curvature and 
displacement effects are present. Additional terms due to the 
variation of the parameters (Ay^ , (fare to be 

included. The multiplicative constants of the eigenftmctions 
are determined by joining the presently obtained asymptotic 
series (for or'-*' ~) with the Euler ized Gbrtler series for the 
blunted wedge obtained in Chapter 4 Thus, the asymptotic solu- 
tion for cr-* provides a test for the assessment of the 
Eulerized series solutions obtained for the nose region, Eurther, 
these two foms can be used to estimate the accuracy of 
numerical methods having uncertain error bounds. 


The material presented in the seciuel is divided into 
two parts. In the first part (6§5. 2-5.4) the general initial 
value problem is treated, while in the second part (6§5,5-5,7) 
we have studied the detailed solutions for a blunted wedge. 

PART I: ‘ Second-order initial value problem 

We consider flow past a two-dimensional or axisymmetric 
body shown in Figure 2.1. The governing equations for the 
momentum transfer problem to be studied here are given in terms 
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of the modified stream functions f (cr , . and l^(cr,r]) by (5.4) 
and ( 3 . 8 ). For completeness, we shall first consider briefly 
the first-order problem. 

5.2 First-Order Problem : 

The first-order houndary-layer equation is 

f ff " +Au^(l-f'^) - 2a-(f'f^- f^jC") =0, 

f(0) = 0 = f ’(0), f'(oo) = 1. (5,1) 

We consider a flow withAy^ = constant = (say) and with 
initial condition prescribed as 

f(<r. ,7?) = (5.2) 

such that f (Tj) deviates slightly from the similar solution 
corresponding to Falkner-Skan solution for p^, known at 
cr- Following Chen and Libby ( 1968 ), the solution for 
f(cr»^) is assumed of the form 

f((r,?)) - fJV) + + f^{cr,r)) + , (5,3) 

where f^ >> >> and so on. Substituting ( 5 # 3 ) in ( 5 . 1 ) 

and properly ordering various terms, we obtain the equations 
for successive approximations. Obviously, the equation for 
foo(^) is the well-known Falkner-Skan equation, viz: 

!H If ^ f o' 

ioo + ioofoo + ^00 ) = 0 , 

fjo) = 0 = fl(0), fl(~) = 1. ( 5 . 4 ) 

I 

The auxiliary condition imposed is that approach of f to 
unity as °° must be exponential. The equation for f^ is 
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found to be, 


tl( 


ff 


“i =fi + v; - 2p„ 

flCo-.O) = fJ(tr,o) = f^(<r,~) = 0, 


2cr (f ’f ' 

' OOT" -j 0- 


If 

f-* f 


oo Icf 


)- 0 , 


(5.5) 


with the initial condition 


m = £*(7?) - ijr,). 

The equations for f^, fj etc. are of the form 

Ilf =s 3} 

== = ^n(“’) =0, 

fjjCcr. ,7?) =: 0. 


Here ia a function of f , , fg-.-f^.,. Using the method of 
eeparation of variables to solve the equation (5.5), we write 


t^{cr,r)) rw, s((r) f(7)) 


Introducing^^ as the separation constant, we find that 

S(cr)-cr n/ while is governed by the following eigen. 

Value problem; 


11 , 

®n + + Cln-2R») Hn + »„ = 0, 

Hn(o) = ir^{ 0 ) = h/(») = 0 . 


(5.6) 


In addition, we have to Impose the auxiliary condition that 
K' approaches zero, exponentially as 7) it la this condi- 
tion Which leads to discrete eigenvalues. It is observed that 
for -.1988..., f^ satisfies the equation (5.6) but is 
not its eigenfunction (since the boundary condition \(0) = o 
is not satisfied). Therefore, we reduce the order oAwo 
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equation (5.6) to two by introducing * Thus we 

obtain 

(pl^)' t (i---tAn)Hjj =0, (5.7) 

where 

p = fj exp ( / dT)) , 

0 

■fc = 

i = i j exp ( / iff) [p„ (3- ^ _ 

The boundary conditions are 

Hh( 0) - [pv£(0)] H^(0) = H„{~) = 0, 

with approaching zero exponentially as 7} -* oo. Equation (5.7) 
is in Sturm-Liouville form and differs from the standard form 
in some details. Using (5.7), it may be shown that for A 0> 
the eigenvalues of the equation (5.7) are real and positive 
(cf. Chen and Libby, 1968). further, it is seen that the ^PT) 
are orthogonal in the sense 

oo 

5^C„. (5.8) 

where C^ is the square of the norm. Similarly, for the so called 
upper branch solutions for -.1988... <Poo < 0, realness ofA^ 
and orthogonality of functions is established directly while 
the fact that the minimum (denoted byX-j) is indeed positive 
is established after it is numerically found to be so. Eor the 
lower branch solutions, equation (5.7) becomes a non-standard 
Sturn-liouville problem and a careful study indicates existence 
of infinite set of both positive and nafffltivp 
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numerical values of the eigenvalues may he obtained by a qbasi- 
linearization technique proposed by Mbby and Chen (1968). 

In view of the above discussion, we conclude that the 
solution for f^ expressed as 

fiCo-,??) = 2 w (7)) (5^5) 

' n=1 ^ 

is meaningful only for the upper branch solutions for 

> -.1988... when positive. Here the is arbitrary 

constant and is to be determined from the initial profile by 
using the relation (5.8), leading to the following expression 
for 

ft (Nyo [(f -fj/fj 47? . (5.9a) 

Substitution of (5.9) in (5.3) yields the solution of (5.1) 
in the following form; 

1(0-, r)) = IJV) + 

+ AgN2(T?)<5-"'^2/^ + ...+ + Ig + .... (5.10) 

5.3 Second-Order Problem : 

Ihe second-order boundary-layer problem is subdivided to 
study separately the effects of the longitudinal curvature, 
transverse curvatiire, displacement and external vorticity. The 
initial value problem for the second-order momentum transfer 
may be formulated in many ways depending on how the initial 
p]cofile is prescribed. One special way is to specify the initial 
profile at a streamwise station (5^ in such a manner that the 
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first-order solution is similar, i,e. f(ot ,7?) = » while 

the second-order profile P(cr^ ,7?) deviates slightly from the 
limiting similar solution However, in general we let 

both the first-order and the second-order. initial profiles 
deviate slightly from the similar profile. We shall formulate 
the problem for the latter case and study the approach of the 
solution to limiting similar solution as (T We describe 
the problem separately for each of the second-order effects. 

5* 3. 1 longitudinal curvature : 


The second-order boundary-layer equation for the longi- 
tudinal curvature problem, in terms of the modified stream 
function , is 





»*• 


(« + 2<r A„ )]/(i ^.A ), 

^ *1 

Cl +Ab^) J{o-, 0) + 2 o-j^(tr,o) =0, 

j'C'T.O) = 0, s'Ca-,?/) „ as 

Here we have defined the operator L''' ae 


(!,£) = I 


III It 

+ f Y 



(2Au^ + E) f Y |.'’y 

^ v + v")- 


(5. 12) 


48 in the first-order problem, we assnme Ajj to be constant 

(denoted byAp alongwlth A„ = The initial condition 
imposed is • 
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such that is slightly different from • In the 

subsequent matter, superscript 1 of f is dropped for conve- 
nience of writing. 

The solution of the equation (5.11) is assumed to be 
of the form, 

■'-SJV) + .... (5.13) 

where << in accordance with the imposed initial condi- 

tion (5.12). Substitution of (5.13) in (5.11) and ordering 
terms, we obtain the equations for ^tc. Alternately, the 

form of F may be written by inspection. Here we adopt the 
latter approach to write, 

F((r, 1?)*^?^*^) + z i. M ■" ^n/2 2 B P (^)cr" 

n=1 ^ “ n=1 ^ 

* ^ ^ Vm ^ ( 5. 14) 

Here Foo(‘^ is the similar solution for specified andA]^. Ihe 
first correction term M^(t 7) arises due to the nonsimilarity of 
the first-order solution H^('^) itself. As mentioned in ^5.2, 
there exists a discrete set of eigenvalues A ^ for the first- 
order problem and this introduces the multiplying factor of 
6' for M^C^) . The next correction term is due to 

the eigenfunctions of the second-order longitudinal curvature 
problem, multiplied by the factors where are the 

corresponding eigenvalues. Other correction terms arise due 
to product of the eigenfunction terms of the f in.at-n-r’H.a-.. ..r.?+v 
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powers and Then follow terns like 

^^ 07 ) with factors which are powers of cT resulting from the 
multiplication of the previously mentioned series, i.e. 

QP- ~(Xn’iym)/2 gQ indefinite proliferation. It may 

he noted that the terms in (5«14) have not been ordered accord- 
ing to the magnitude of the power of cT j such ordering is done 
after specifying the values of the pressure gradient parameter 
and longitudinal ciu-vature parameter The multiplicative 
constants are to be determined from the upstream data for 
the longitudinal curvattire problem while the constants are 
known from the first-order problem. 

Substitution of (5.14) in (5.11) and collection of 
coefficients of the like powers of cT yields the follovang 
sequence of equations: 

C+ 

+ (2poo + > 

I Jo) = 0 =» I'l.(O), f 1(7?) ^ -7] as ry- oo. ( 5 . 15 ) 

I<(Pn» Ai-Vj^) =0, (5.16) 

P^(0) = P;(0) = P;(oo) . 0. 

I-(M„,Ai-Au) = ‘-(^■»Ai)PooN" t (2 Poo+Ai-A„) 

+( i+PJn^"+(A3.-Poo- 1) CVfX 

+ (2R=+Ap 
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+ (2p^ +A3_-2An)-^n ^'r 

Vn [nX -X,Vn 


■4An^n » 


, ^ ''rpn-'n 

Qjj(o) s= QjjCo) = QjjC"’) = 0. 


Here I 1 b the differential operator defined as 


(5.18) 


^(^! S)_ Y + f _ ( 2 p^ + t) f y + (l+e) f^'y ( 5 . 19 ) 


The first equation (5.15) is the familiar equation which yields 
self-Blmllar solutions. Ihe next equation (5.16) is a homogene- 
ous equation with homogeneous boundary conditions, leading to 
an eigenvalue problem. If we put y^. ^X* in equation (5.16), 

4 4-- Ti I 


it reduces to 


^n'-^ Vn - (28.- A*) fl i. (1-A*) f^' p = 0 


’^n(°^ = ®n(°> = ^n(“) = 0- 


oo 'xi 


(5.16a) 


The form of the above equation is identical with the equation 

( 5 . 6 ) for H^(1?) whose solution is already known. The rest of 

the equations for M^(7?) , q^( 7)) etc. are to be solved by numerl- 
cal integration. 


5 . 5.2 Transver se curvature: 


is given by 


The governing equation of the transverse curvature problem 

..... T- 


I* (•I’^'*'! Ag ) = - (2L + f'") + f" + ff' ^ 

• * 4". c 


(1 +Ab ) F^'^)(cr,0) + 2cr'F^)(cr, 0) = 0, 


w 

j'('‘\o-, 0) = 0, J''(*\o-,17).-7? as7?-». 


(5.20') 
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Initial condition: 

The transverse curvature parameter is assumed to be constant 
and denoted by alongwith Ay = Poliowing the arguments 
used earlier for the longitudinal curvature problem, we assume 
the following expansion for 

P((r,??)- PJT?) + 2 A M (7))0''-‘^n/2 ^ ^ B P (7?)cr -^n/2 

n=l ^ n=1 ^ ^ 

+ 2 cr^^4-... (5»21) 

n=1 ” 

Here we have dropped the superscript t from F. Substituting 
(5.21) in (5.20) and collecting various powers of cr' , the 
resulting equations are written- in the operator form as given 
belov/, 



. ft? 

A^) = -T)(2p„ + fco ) + i 

•OO + foo^c» 


S' Jo) 

= 0=jI,(0), rlra-7) 

as 7)-* °o. 

(5.22) 


o 

If 

1 

+5 

< 



P„(0) 

= p;(0) = ?;<») = 0. 


(5.25) 

KMn. 

=-(i+VpX'-*- 

(2 Poo 



-d-v 

♦It n 

vfLK 


+(i-V 



Hjj(O) 

= M^(0) = = 0. . 


(5.24) 


I'(Qn> At^^Ajj) (l+A^~Ajj) ^xi^n ^ 
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Equation (5.22) yields the limiting similar solution. With the 
substitution of ~-^n’ equation (5.23) reduces to 

(5.6) whose solutions are known. Other equations are to he 
solved by numerical integration. 

5.3.3 Displacement speed ; 

The second-order boundary-layer equation due to dis- 
placement is written as 

L-** (E^^>,Ag ) = -Ab - 2l\„ , 
d "d ^1 

(I^Ab ) (o", 0) s- 2cr F^^(<r, 0) =0, 

d ^ 

F’^'^^((r, 0) = 0, <->):= 1, (5.26) 

Initial conditioh; 

P^‘^^(ot,??) = ^""(1?) 

Here again A b^ = constant =A^(say) andAy = As in the 
case of the longitudinal curvature, we assimie the following 
expansion for after superscript d has been dropped 

for convenience, thus 

F((r,-3?)-' EJT)) + s AM (7?) cr ^ ^ B P (7?)cr 

n=1 ^ ” n=1 “ 

^ + ••• . (5.27) 

Substituting the above expansion in equation (5.26), we obtain 
the following set of equations in the operator form; 

FjO) = 0 = f 1(0) , Fjoo) = 1, 


(5.28) 
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= 0 


_ Pj^(O) = \(°°) = 0. 

I'(Mn»Ad“V = -(''+A^)^'co \ +(2p^ 

M^(0) = M^(0) = M^(~) = 0. 


(5.29) 


(5.30) 


I'(Qn» Ad-2An) = ~(''+Ai-An^ K\ 

-( 1-\) Vn ’ 

%( 0 ) = Q^(0) = Q^(oo) = 0. 




(5.31) 


Equation (5.27) yields the limiting stailar solution for the 
displacement effect. Again we can substitute 7^^ +A in 

the eciuation (5.29) to reduce it to the form (5.6) for which 
the solutions are known. The other equations for etc. 

are integrated nmerically. However, for the special else of 
Ad = 0, the equation (5.30) for M^(7?) has an exact solution in 
terms of the first-order eigenfunctions N^(^) given by 

Mn(^) = [<(^) ■^hy^]/2. • (5^32) 

5.3.4 External vorticity ! 


The second-order boundary-layer equation for the external 
vorticity is given by 

1-2V) = -a - So-oc^, 

2 (l-Ao ) I't''^(<5-,0) + 2<3-p('") (cr, 0 ) = 0 , 

1 ^ 

E'^^^(cr, 0) = 0 , ?'^^)((r,7?)- as 7}-^ 00, 


(5.33) 
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Initial condition: 

Here again, we assume = Poo = constant and drop the super- 
script V from F, As for the otlfer second-order effects, here 
again we assume the following form for I': 

I'((r,7))^pj7y) ^ 2 A^M^('>?)a-^n/2 ^ 2 B P (T/) cr- -Xi/2 

. (5.34) 

Substitution of (5.34) in (5.33) yields the following system of 
equations :• 


I'C^’oo. 1~2pJ = - liifl (7}-fJ7j)) 

'j^oo 

j’jo) = 0 = pjo), I'jrD^ri as ??- oo. 

“-Vn) = 0, 

^n^°) = = 0. 

I-CMn, 1-2p^ = -(2-2pJ 

Mj^(O) = M^(0) = M^(oo) ^ 0. 


(5.35) 


(5.36) 


(5.37) 


L(Qjj, 1-2poo -2An) - -(2--2p^-> ) N M +(1-2A )n' m’ 

11 XI n n' n n 

-( 1-;^) <’ , 
v°> = e^(o) = q'(=o) = 0. 


(5.38) 


The limiting similar solution due to -external vorticity is 
given by equation (5.31). Substituting Y _i4.2p^ =A* in 
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5.4 lumerloal Results and Discussion : 

Ror each of the second-order problems along with the 
first-order problem, detailed solutions have been obtained for 
two representative values of namely 0 and 1/2. The effect 
of variation of the second-order parameters, arising due to 
longitudinal curvature, tra,nsverse curvature and displacement 
speed on the corresponding solutions is examined for 
-2 <4. For completeness, we first discuss briefly 

the first-order solutions for = 0 and 1/2. 

The solution of the first-order problem can be expressed as, 

f(0-,7?)- fJ7?) + A,It^(7?) or-V2 + A2H2(7J)cr-V2+... .(5.39) 

The eigenvalues and eigenfunctions are determined by 

solving equation (5.6) with the auxiliary condition specified 
earlier. Libby and Fox (1965) have listed the first ten eigen- 
values for the above equation for = 0 while Chen and Libby (1968) 
have provided the first twenty eigenvalues for = 1/2. The 
important question that arises is where to truncate the series 
(5.39) for the purpose of studying its effect on the second-order 
solutions, . This question is linked with the prescribed initial 
condition; for a very small departure of the initial profile 
from the similar one, only a few terms may provide a satisfactory 
answer. Since, we are interested here in small deviations of 
the initial profile so as to linearize the solutions about the 
limiting similar solution, it is decided to include terms upto 
cr“''^2/2, This will be consistently followed in the second-order 
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series representation for I'(cr,')7) also, The required, first two 
eigenvalues and the corresponding squares of norms for = 0 
and 1/2 are given in Table 5.1. 

Table 5.1i The eigenvalues and squares of norms for 
= 0 and 1/2. 



n 

n 

°n 


1 

2.0 

2.2672 

0 


2 

3.774 

1 . 1447 


1 

3.092 

2.066 

1/2 


2 

4.959 

1.011 

The multiplicative constants A^ and 

A 2 are 

estimated from 

relation (5.9a). The 

eigenfunctions 


are shown in Figure 5.1. 


It may be noted that all the eigenvalues are positive so that 
the initial profile which deviating slightly from the similarity 
profile, approaches asymptotically the similarity profile as 
cr-* Vie nov7 proceed to examine the corresponding second-order 
boundary layer results and see if such a behaviour is exhibited, ■ 

The solution for the second-order problem is expressed as 

+ E + ... ^ 

n 

(i = l,t,d,v) 


(5.40) 
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Two important points to be considered in regard to the above 
solution are: to choose a criterion for truncating the series 
and to order the terms according to the povrers of <r . For the 
first question, we mentioned earlier while discussing the 
first-order problem, that we shall retain terms upto cr’”‘^2/2 
only so as to make f(cr,'?7) and F(cr,')7) consistent in some sense. 
Thus, in the follov/ing discussion we shall show only those 
terms which have the magnitude of power of CT less than A. 2 / 2 . 

Kov/ to answer the second question, we remark that, in general, 
it is not possible to order the terns without specifying the 
parameters PogjAijA-^ andA^. This is due to the dependence 
of Poo and of A^, A^ and We now discuss each 

of the second-order effect separately to show some typical 
properly ordered solutions for F^^^(o^,77) for various values 
of the parameters, 

5 . 4.1 Longitudinal curvature : 

The longitudinal curvature problem is governed by the 
equations ( 5. 15)-( 5. 18) . The eigenvalue problem (5.16) involves 
the eigenvalues which is determined by using the relation 

y„ = (5.41) 

and in particular 

= Ai +A., , Yg = +A 2 . (5.41a,b^ 

Depending on the longitudinal curvature of the body and hence, 
the value of A the terms in (5.40) can be ordered according 
to the powers of cT ; five such distinct cases (a-e) given below 
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(a) Ai > U 2 “Xi> 

n<r,rj) ^ fJt?) + A^M^(r/)cr-V2 ^ j,^^^r))c--^2/2 

+ 0(cr-^2/2) 

(b) (Ag-A^) >Ai > 0 

.+ A^M^(77)(r'-V2 ^ B^P^(>?)cr~V2 

+ A^l^m cr'^2/2^ 0(6" "^2/ ^) , ( 5 . 45 ) 

(c) 0 > Ai > -(^ 2 -Xf) 

i'’(^,??)-FjT9) + B^P^(r/)cr-'^l/2 A^M^(i7)(r"*V2 

+ BgPgC??) cr "*^ 2/2 ^ _>2/2 ^ o((r ''^2/2) 

(5.44) 

(d) -(A 2 -A 1 ) >Ai > -A^ 

F(<r,7?)-Pj7?) ^ B^p^(7J)cr“%/2 ^ B 2 P 2 (^/)cr"V 2 

+ A^M^('r))cr “^ 1/2 + A^2('^)cr 'V2 ^ o(^-.A2/2) 

(5.45) 

(e) Ai < ~Ai 

P(cr,7])^B^P^(7?)(r-~V2 ^ bjtT) + B^^{r})^-y2/2 . 

+ A^M^(7))(r "^1/2 + ^^^(^)^-A2/2 ^ q(^-.;^2/2) 

(5.46) 

It may be noted that we have intentionally omitted the 
equality .sign from all the above shdwn limits for ■^]^5 the 
reasons for doing so will become apparent later where we 
discuss the solutions corresponding to the values of A^ 
given by equality sign in the above mentioned limits. 
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We first discuss the last case (e), since it has certain 
implications about the nature of the longitudinal curvature 
of the body. For this case, the first eigenvalue A ^ is negative 
so that the corresponding term due to the first eigenfunction 
P^(7?) is multiplied by positive power of cr . Thus, for 
Af < A-] > "the solution is unstable in the sense that any depar- 
ture from similarity (at a given initial station) will grow 
indefinitely in the downstream direction as — o®. If, however, 
•the initial condition is prescribed such that = 0, then 
the equation (5.46) shows that the solution is again stable, 
in the sense defined above, till A is decreased so much that 
Ai < - A 2 so that ^2 becomes negative and again makes the 
solution unstable. Similar arguments may be applied to show 
that for various prescribed initial conditions, - J\^ 
etc. are the critical values. For a general initial condition, 
however, it is the value A^ = itself that is the most 
critical and solutions will be unstable forA^^ < -A^. We 
recall here an earlier discussion of §3.3; it was shown for 
the longitudinal curvature problem that the self-similar solu- 
tion, F^('^) in the present notations, has several infinite 
discontinuities whenever A -j_ takes a value equal to the eigen- 
value A^^fo^c- given p (poo present notation). That the 
locally similar solutions are indeed singular at these criti- 
cal values of A]_!i is further reinforced from the present 
analysis which shov/s that the flow becomes unstable for the 
values of Aq^ corresponding to the singular points. In §3.3, 
while discussing the design of aerodynamic shapes, it was 
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suggested that the curvatiu?e of the body is to be so pres- 
cribed that Within the linear theory, the present 

♦ 

analysis, puts a further restriction on the selection, of the 
value for I't seems desirable that the values of A 2^5 such 

that ^ » should be avoided to ensure that any slight 

departure from the similarity profile will not cause the skin 
friction to grow 'indefinitely in the downstream ■ direction. 

Por the values ofA^^ shown in other cases (a~d) , the 
solution (5.40) is stable in the sense that the slightly non- 
similar initial profile does approach the limiting similar 
profile as o ' The governing equations ( 5. 15)-(5» 17) are 
integrated numerically for various values ofA^ and Results 
are presented for two representative values of = 0 and 1/2 
and A 3_ in the range -2 < Aj^ < 4. The solution of the equa- 
tion (5. 16) and (5.17) is obtained for n = 1 and 2. Figure 5.1 
shows the eigenfunctions P.j('^/) and for both = 0 and 

1/2. The corresponding eigenvalues y and Yg for a given 
are given by the equation (5.46). 

The solutions of the first two equations of the system 
(5,17) are shown in Figures 5*2 and 5.3 for tw/o values of 
.i.e. 0 and 1/2. Figure 5.2 shows M^'^^^(O) vs. A^^ obtained 
from solution of the equation 5.17 for n = 1. It may be seen 
that several discontinuities occur for both the values of p^ 
at certain critical values of the longitudinal curvature 
parameter A^^. First discontinuity occurs at A 2_ ~ ihis 
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behaviour of the equqtion (5*17) is easily explained by consi- 
dering the operator of the equation v/hich, on Substitution 
of =0, reduces to the first-order eigenvalue operator 
(5.6) corresponds to Thus, becomes the eigenfunction 

of the operator of the equation (5.17) and leads to the dis- 
continuity (nonunique solution) in the numerical integration 
of the equation. The next discontinuity occurs at A ^ =A^“A.2’ 
Again, substitution of the value reduces the operator of (5.17) 
to the first-order eigenvalue operator (5.6) corresponding to 
Xg snd hence, the singular behaviour. There will be many 
more such discontinuities whenever A^ fakes a value equal 
toXi~Ajj» A “I* solutions for are also found to 

be singular whenever A]_ = -Ayi* This happens due to the 
singular behavio\u? of the corresponding equation (5.15) as 
was explained in Chapter 3 (§3.3). Thus, for these values, 
the Right hand side of equation (5.17), involving .the solution 
of (5.15), blows up. 

The solution of the equation (5.17) for n = 2 is shown 
in Rigure 5.3 in the form of M2^^^(0) versus A Here again, 
several discontinuities are observed whenever A^^ ===A2-An 
ox A-^ = “An reasons mentioned in the preceding 

paragraph. We may note that the first singularity occurs here 
for positive values of which has never been observed before. 
Therefore, there exist both positive and negative critical 
values of Aj^ sT^ch that the second-order longitudinal curvature 
solution becomes unstable in the sense that any departure from 
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similarity at a given streamwise station grows in the increasing * 
streamwise distance. 

As mentioned earlier, the equality sign was purposely 
omitted while writing the limits on A i for different cases 
(a-e). It is now clear that for those limiting values of A ^ 
with equality sign, the flow is unstable because the solutions 
and are themselves found to be singular at these values, 
lo what we said earlier about choosing the values of A^ for 
aerodynamic shapes, it may be added that besides the restriction 
>■' the isolated critical values of A-j^ (given in the 
preceding paragraphs) are also to be excluded, 

5 . 4 • 2 Transverse cin :va ture; 

The . transverse curvature problem is, governed by 
equations ( 5. 22)-(5. 25) * The value of is detemined by 
using the relation, 

Yn " ^t ■•‘^n "" 1,2,...) (5.67) 

As in the case of the longitudinal curvature, we may order 
the terms in (5.40) depending on the value of A^ and obtain 
the five distinct cases given by equations ( 5 . 42)-( 5 . 46) after 
replacing Ai byA^* ^’or any givenA^J nature of the solu- 
tions being similar to that for the longitudinal curvature 
described in §5.5.1, no further comments are needed. 

The solutions of equation (5.24) for n = 1 and 2 are 
presented in Figures 5.4 and 5.5, The variation of M^^’'^^(0) 
and M 2 ^'^^( 0 ) with for both = 0 and 1/2 is shown, 
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zespecljively y in Figures a onri k c < 

feures 13.4 and 5.5. Again, several singula- 
rities are observed for +Via t 

critical value^ which are 

the critical r\-F A • 

cai values of given in §5.5.1 and, 

therefore, no further comment is needed. 


5.4*3 £i L.A I 

Solutions to the displacement speed problem are obtained 
by aiumerioal integration of the equations ( 5. 28)-(5, 50) , The 
Value of is again obtained by using the relation 


Xi ■^d '’'‘^n “ I52 ,..,), (5,48) 

The terms in the expansion of are ordered 

according to the value of as given by the equations (5.42)- 
(5.46 ). The comments given earlier regarding the ’ stability of 
the flow for various values of the second-order parameters 

AijA-fc hold good for the displacement effect also and need no 
elaboration. 


The solution of the first equation (n=1) of the system 
<5.30) is shown in Figure 5.6 for two values of namely = 0 
and 1/2. Again, we observe from Figure 5.6, that m"(o) has several 
singularities forA^^ ~'^1“An (^>1)* These infinite discontinui- 
ties occur at the same values of the displacement speed parameter 
A^ as observed for the longitudinal and transverse curvature 
problems but for one exception: at A^ = 0, the solution is finite 
and is marked by **' in the Figure 5.6, For the last case 
(A^ = 0), the operator of the equation (5.30) is reduced to the 
first-order eigenvalue operator , suggesting thereby that 
inspite of (at A ^ = O) being an eigenvalue, the function 
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is not singular. The particular integral of (5.50), viz. 

satisfies all the houndary conditions, as does the complimen- 
tary function 

Consequently, v/e obtain the following expression for the com- 
plete solution of M^(17): 

m^iV) =A.M**(7)) (5.49) 

where A is an arbitrary constant. One v/ay to evaluate the 

4. 

constant is to join smoothly the two solutions for A ~ ^ 
and 0“. 

The solutions for function 12(0) as a function of 
are shown in figure 5*7 for the values of 0 and 1/2 » Again 
we observe that the solutions of the equations (5.30) for. 
n = 2 are not unique whenever ~^2'"^n ^ S'igure 

5.7, forA^ = 0, we have again indicated the corresponding 
value of 12(0) by The nature of the solutions is similar to 

that of and hence no further comments are needed. 

5.4.4 Internal vorticity ; 

The second-order external vorticity problem does not 
contain any new parameter. Solutions to equations (5.28)-(5.30) 
are obtained for - 0 and 1/2 by numerical integration. The 
V3,lu0 of is olDiainod fey ifeo rolaiion 

= 1-2p„+An (n = 


( 5 . 50 ) 
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Ihe terms i„ the expansion of oan be ordered for a 

given It May be noted that for = 1/3, equation (5.30) 
ie reduced to the foi» same as that for the first-order 
eigenvalue problem so that the solution of M^(ri) becomes 
singular for all values of n so that the flow is unstable in 
the sense defined earlier. For = 0, the terms in (5.27) 


are correctly ordered in the following manner; 

C. PJT?) + i^m,(7)) cr-^1/2 + B^P^(0j)cr-V2 

+ A2Mg(l))cr-V2 + 0(cr-^2/2)+,..(5.5i) 
The operator of equation (5.30) for n = 1 and = 0 (with 
corresponding X .j = 2) is reduced to the following form; 

.1)^ m’” + fjs” + 


This leads to solution M^'(O) = 0 in the numerical integration 
of equation (5.30) for |3^ = 0 and n =1. The solution of 
(5.30) for = 0 and n = 2 is found to he (O) = -2.591. 

Substituting the values of M., and thus obtained in (5.51), 
it may seen how approaches to the limiting similar solution 


for p^.= 0. 


Part II; Skjn friction on a blunted wedge for large 
streamwise distanc es ; 


The general analysis of Part I essentially deals with 
the evaluation of the eigenvalues and the corresponding eigen- 
functions which are to be included in the expansion of the 
Variable that represents the solution to an initial value 
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problem of the type studied in the preceding sections. However, 
m many problems of praotioal interest, the parameters A„ 
etc. are,ln general, functions of o-and, therefore, the e^pan-"- 
ion must include additional terms (besides eigenfunctions) 


now 


to account for the variation of A , A etc. with o- . We 
obtain detailed solutions for one such problem i uniform flow 
past a blunted wedge shown in Figure 4.1, and study the momentum 
transfer problem for large <r. . The corresponding solutions for 
™all values of were studied in 44.4 of Chapter 4 and the 
Eulerized form of these solutions is used to impose the initial 
condition for the present problem. Thus, we shall show how the 
Eulerlzed series solution for nose region approaches the 
limiting similar solution at o'- quantity of particular 

interest is the skin friction coefficient defined by (4.55) and 
as such we focus our attention on f "((r„0) and p'W, 0) defined 
(4*54) and (4*55) respectively, 

5 • 5 girs t-Order Prt^-hlpni ; 

The first-order boundary-layer equation in terms of the 
modified stream function f(<r, 7 )) is given by (5.I). The initial 
condition (5.2) is prescribed in the following form: 

f ((T^ , 0) = (5.50) 

A 

where defined by (4.85), is the Eulerized Gdrtler series 
solution obtained in Chapter 4. The pressure gradient parameter 
Au for the blunted wedge given by (4.21), viz. 


Au = (l“Pco + cr)/(l-p^ + 2cr), 
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is expanded for large values of cT to yield the following expan- 
sion: 

Au. = Poo + S (Poo-1)”^^ (2cr)“^ . (5.51) 

1 n=1 

Using (5.10), the solution of (5.1) is expressed in the follow- 
ing form: 

= f J7?) cr"'' + 

+ + A2N2('?)(r-V2 + ... ,(5.52) 

Here the additional terms of the type are 'forced' 

by the pressure gradient parameter An since it varies with cT 

1 ’ 

in the manner given hy (5.51). It is noted that the relative 
position of the terms in (5.52) will depend on the value of 
In (5.52), v/e have shown correct order of the terms for 

0 < Poo < 1-. 

Substitution of (5.52) in (5.1) and collection of the 
coefficients of the like powers of cr yields the following 
system of equations (bar on f^'s is suppressed for convenience) 

f.j + f^f., + 2(1-PJ foof^ - = -d-poo) ( 1-^00 )/2» 

f^(o) = f^(o) = f](°°) = 0. (5.53) 

^2 ^ ^00^2 ^ 2(1-Poo) ^00^2 "" ^2 ~ ^1^1 2-Poo) 

+ d-pj^ V d-PJ^d-f^A^ 
f^io) = 4^0) = f2(“) = 0. (5.54) 

Ihe equations for fj"'?) and N^(7?) are same as those given by 
(5.4) and (5.7) respectively. Ihe solution of (5.53), expressed 
in terms of the solution of the equation (5*4), is of the 


following form: 
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= (1-Pj^ (40"'' (.VL- fj. . (5.55) 

The solution of (5.54) has to be obtained by numerical 
integration. The multiplicative constants Ai , A 2 etc. in 
(5.52) have to be determined by making use of the initial 
condition ( 5 . 50 ) whereby we join (5.52) to the Gdrtler 
series solution obtained in Chapter 4 . 

5.6 _Second--Qrder Prnbipm i 

The second-order boundary- layer equation for the 
longitudinal curvature is given by (5.11) and is repeated 
here (superscript 1 is dropped) for convenience, 

Ab ) = [(Ab +An - 1)(f" + ff' + 2(rf'f_ 

“ ! 

4o” j f f d^/) + (2 Atj ’^'A'r ) 

7] 1 

(a + 2(r + A u^)]/(1 + Ab^) , 

(l+.'\Bi)^'(cr,0).+2crI’^.(cr, 0) -0 = ^'(cr, o) ,I'’(o- 7), ( 5 . 56 ) 

as 

The operator l'*’ is defined by (5.12). We may note that A ^ 

# 1 
and A-r are function of 6~ here. The initial condition is 

■“1 

prescribed by fixing the initial profile at any given station 
cr^ in terms of the Bulerized Gbrtler series for the longitu- 
dinal curvature, obtained earlier in Chapter 4 . Since we 
are primarily interested in the skin friction coefficient, the 
initial condition has the following form: 

B* , 0) = (5.56a) 

whore is given by (4.87). The quantity defined by 

( 4 . 56 ), is expanded for large values of o-'to yield the 
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following expression for denoted by B*: 

\ = (2cr)-^ ( 2 o -)~2 ^ (5.57a) 

Using ( 4 . 37 ), we obtain the following expression for the 
longitudinal curvature parameter for large values of cr : 

Ag = -2 ^ Z (-2) (2(r)“^ =A* 

1 n=1 J-o 

+ ... . say (5.57b) 

The above form of the expansion suggests that the expansion 
of i!^ \o will contain a series of functions multi- 

plied by inverse integral powers of cr . Including such 
additional teims in the expansion ( 5 . 14) of I'(o",'''n, we 
write the solution of (5.56) in the following form: 

F(<r , 17 ) - fJ'' 7) + B^P^(T/)cr "'^' 1/2 ^ f^(r/)cr''+ 

+ A^M^('^)o- ■"■'■1/2 + 

. (5.58) 

1 2 

Here tenns in cr and cT arise due to variation of A -d 

Bf 

with o~ . The relative position of the terms in (5.58) will 
depend on the value of in (5.58) the terms have been 
ordered correctly for = 1/2. The value of is determined 
from the relation (5.41). 

Substituting (5.58) in (5.56) and collecting the like 
powers of o'", we obtain the following system of eq.uations for 
the longitudinal curvature problem (again, bar on P^^'s is 
suppressed) j 
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pj= v: .( 2 ^ 


+ ( 1+A f ) f " S' _ 1 

± ^ oo “*-00 “■ 

0 


IfT 

-Vi„ + 2(P„-I)(p„+1)-1 (c(„+i7pj 

+ (P«-3)(Pco+l)’''(fJl'+f^l), 

Jjo) = 0 = FJO), F (17) ^ (5.59) 

^2(^1- Ay pj = j;"+ f^;' (2_2 p^ .y^*^) 

+ (Ay1) fl'F, - , _° . y 

4-0 T oo 


‘n **' -1 

-'?f^ + 2 (g„- 1 )(p^+ 1)-1 [-“-l+^C P„- 1 ) ^ 2 ] 

+ (p„-3)(p„+i)-^ [f; * 

+ '/fl^/2 - (l+pj'’'(fr+f„fl + ao+r)pJ/2j 

+ (P«r3)CP<xr'l)(Poo+'*) ^(fco + ^«^00 + ttQ +l)Poo), 

F,(0) =l-;(0) =?'(«.) =0. (5.60) 

S^“i’Ab,’ P«) - “i + + (A^-sp^-A* ) fjil 

X Q I 

+ (1+AyAi)f»M^ + (2p„+AyAi) 

H>l+(1-Ai) = -’’<'-2(A»-i)(Pc„+i)-^(i-Ai) 


[^i^^.fi(i7)] * (p„-3)(p„.n)-1 p; + (1-Ai)fjj, 
+ + Ai{i+P„-Ai)“’'i + fjjJ + (1-Apfjj^ 

+ (l-Al) ( llm M,(7?)]] , 

' rj-^ oo * 

M.j(0) = m !^( 0 ) = M.j(°°) = 0. 


(5.61) 
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PJ3 p;\ * (y^-2p„.A* ) fV 

o 

+ <1+A*^- y^) = 0, 

i’i(o) = p!^{o) = pj(~). 

Here 


(5.62) 


a 


'o 


= lim [r,_ 

' / -♦ oo 


“1 = [ -f1«)], 

fj^oo ' 

■sf 


a: = -2. 


A 


0 

* 

1. = 


^oo* 


The second-order boxmdary-layer equation for the dis- 
placement effect is given hy (5.26) and is repeated here 
(superscript d is dropped) for convenience, 

1-^(H,A3^) = -Ab, 

(5.65) 


d ^d. 

(1 + Ag ) i'((r , 0) + 20 - cr , 0 ) = o 


^'(o', 0) = 0, F'(cr,7)) = 1 as "y; - oo. 


(5.64) 


The initial condition; 

, 0 ) = 

Here Tg ^ is given hy (4.79) and represents the Eulerized 
series solution of the displacement problem, obtained in 
Chapter 4. The quantity B^ [ = U 2 (s, 0)/U^(s, o)] , defined by 
(4.51), is expanded for large cr to yield the following 
expression, denoted by B^; 

= -(1--pJ(2cr )~'’/2 [Bq + cr-^ + ... J 




155 


where 

Bq = ( 2 a ^ - a .^)/2 

= (1 - U(6a^ ~ 7a^)/8. (5.65> 

Here a^, a^ etc. are same as those defined by (4. 50), 

Using (4.32), we obtain the following expression of the 

displacement speed parameter A -0 for large cr : 

d 

= -1 + (i_g^) [(6a^-7a^)/2(2a,-a„)] cr -I + .... 

^ ^ *1 

^ ‘*‘^d ■** (5.66) 

0 1 

Ihe five term series (4,31) of was obtained starting from 
the stagnation point and has been utilized above to obtain 
its asymptotic form ((T-*°=>) with the hope that the factor 
(.1 +% extracted from the series gives the correct 

behaviour of B^ [and, therefore, of U2(s,o3 as o'- 00, 

In the strict sense, we should find the required expansion 
for large <r , starting from I /0 = 0. However, it is still 

instructive to use the expansion (5.66) and obtain useful 
estimates of the solution of the displacement speed problem 
for large o'*. 

The solution to the displacement speed problem for 
large cr is assumed to be of the following form: 

l'(cr,7?) =l'jm + +B.,p.,(r/) <3--V2 

+ "'' 2/2 + . • ■ . (5. 67) 
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Here the terms of the tvop ? (TS\nr ^ o. 

o Tiype i!^^//;cr arise due to variation 

^ith 0“, given by (5.66). The relative position of 
the various terms in (5.67) are shown correctly for = I/2. 
The value of is to be determined using (5,48) with 
Ad-^A^^. 

Substituting (5.67) in (5.63) and collecting the 
coefficients of the like powers of 6~ , we obtain the following 
system of equations (again, bar on P^'s is dropped): 

a 

*'0,(0) = 0 = F^(0), F„(~) = 1. <5.66) 

Pj = -(l-pJ(6a^-7a„)(2a^-a„)-V2-(p„-i)®. 


^.((O) = F^(0) = P^(oo) = 0. 

(5.69) 

f Poo) = o> 

' d 


M^(0) =m'(0) =m’(-) = 0. 

(5.70) 

Poo) “ 


P.j(0) = p’(o) = p’(°o) = 0. 

(5.71) 


Here 

A^o = -1' 

Aj = (1-pJ(6a.|-7ag)(2a^-ag)~V2- (5.72) 

The operators , Lg, and are same as defined earlier for 
the longitudinal curvature problem. 
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5-7 Em r xc ^l Solj^ ^ and D.i3nn.c,...-inv, > 

In order to examine how the skin friction will 
approach the limiting similar solution value as cr-* oo^ we 
now evaluate the asymptotic expansions (5.52), (5.58) and 
7) .and estimate the multiplicative constant 'of the first 
eigenfunctions. The procedure will he illustrated for one 
typical value of namely = 1/2, 

The first-order problem has been studied by Chen et.al. 
(1969) and numerical solutions are provided for p^ = 1/2; 
for completeness, wc mention briefly their main results. J'or 
Poo = 1/2, the first five terms of (5.52) have been calculated 
to yield the following asymptotic expansion for the first- 
order skin friction: 

f'*((r, 0 )- 0.927680 + 0.115960 0 ““'' + -1.546 

-0.02615 O' “2 ^^^-2.46 (5.75) 

The value of A^ is determined by fitting (5.75) to the 
Eulerized Gbrtler series solution, given by 

= 1.232588 - 0.246920Z - 0.033951z^ -0.01081 6z5 

- 0.0047532^- 0.002510z^. (5.74) 

This leads to A^ = -.026. Chen et.al. (1969) have also used 
an alternate approach whereby they utilize the orthogonal 
properties of the eigenfunctions and find the value of 
A^ = -.027 which is close to the value obtained by the first 
method. Substituting the value of A.| in the series (5.73), we 
can find how the skin- friction approaches the limiting Ealkner- 
Skan value as (y- 
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We now proceed to determine the multiplicative constant 
of the first eigenfimction corresponding to the second-order 
problem. Let us first consider the displacement speed prob- 
em. Equation (5»71) is solved subject to the normalizing 
Gondi-tion P^(o) = 1. g)he first set of eigenvalues is given 

Kn ^^n for = 1/2 and n = 1 ; it leads to the 

value = 2.092. Equations (5*68), (5.69) and (5.70) are 
numerically integrated and we obtain the following solutions.* 

^00 (0) = 0.5389, (0) = -2.7668, E^*(0) = -9.0579. 

Substituting the above values along with P^(o) = 1 in (5.67), 
we obtain the following asymptotic expansion of 0) for 

Poo = 1/2: 

S'’’ (d", 0) V.’ 0.5389 - 9.0579 cr-'^ + + 0(cr • 546j ^ 

(5.75) 

The multiplicative constant can be determined by using 
the initial condition (5.64). However, as the asymptotic 
expansion of B^ starting from 0*“^ = 0 is not available, we 
have used for the asymptotic expansion the same series of 
as in the case of the Gdrtler series for small 0*^. We, 
therefore, use the initial condition in terms of the Bulerized 
Gbrtler series of F'*((r^, 0) equated to asymptotic form of 
S' f 0) given by (5.75). For = 1/2, the Eulerized G6rtler 
series of F (cJ^, 0), obtained in Chapter 4, has the following 
forms 

F"((r, 0) = 1.8488 - 1 . 7322 Z - 0.7194z^ - 0,4978z^ 

- 0. 45672^ - • • • , ( 5 . 76 ) 
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where z = cr [0-+ (l-pJ/2]“^ 

The outline of the procedure for obtaining is as follows: 
yfe rewrite the first three terms of {5*15) as a series in 
( 1-z) , thus 

r'^iO) =-35.6926 + 36.2316 z + [4( 1-z)] ^ 

+ 0[(1-z)■'•5'^S]. (5.77) 

Eq.uating right hand sides of (5*76) and (5*77)j we obtain 
the following equation for B..|i 

B^ [4(1-2)]^*°'^^ = 37.5415 - 37.9638 2 - 0.7194 

-0.4978 2^ - 0.4567. 

= 37.5415 (1-1.011 2 -.0191 
-.0132 2^ - .0122 2^'). 

» 

Raising both sides to the power (I/I.046) and using the 
binomial expansion for the right hand side, we obtain the 

following series for B.^ : 

(3 )1/1.046 4(i_2) = (37 . 5415) (1-.966 z - .0183 

- .0126 z5 - .0115 zb (5.78b) 

Again, the right hard aide of (5.78b) la rewitten In powers 
of (1-z). Thus, 

(3 ,1/1.046 ^ (57. 5415) [,(1-0.9660-0.0183 

-0.0126-0.0115) + (0.966+0.0366 
+0.0380+0.0463) ( 1-z)+0( 1-z) . (5*79) 
As it is, the two sides of the above equation do not match. 



140 


But on the right hand side of (5.79), the first coefficient 
in the square bracket has the last partial sum equal to the 
value -.0095 and as we would expect, it seems to approach 
zero. Therefore, we put the first coefficient equal to 
zero so that the coefficient of (l-z) can be matched. The 
second coefficient has the last partial sum equal to 1.087. 
Thus, equating the coefficients of (l-z) on the two sides of 
the equation (5.79) we obtain, 

(57.5415) ( 1. 087/4) 

= 9.57 

Substituting the value of B^ in (5.75), we obtain the following 
series for B ^‘^^(cr,0): 

b"^^^(o-, 0) .5589 - 9.059 cr + 9.575 

+ 0 (cr“'' • 5^^^) . (5.80) 

To obtain the asymptotic series for the overall displacement 
speed contribution to the skin friction, we multiply (5.80) 
by the asymptotic form of B^^ (5.65) to yield the following 
series: 

B* f , 0) - -0.017 + 0.289 cr 

+ 0(o-'~'^*^'''^). (5.81) 

The above asymptotic series for the overall contribution of 
the displacement speed to the skin friction shows that as 
cr the contribution is zero; thus it confirms what was 

predicted earlier (see Figure 4.11) by the Bulerized Gortler 
series solution and, in addition, gives the rate of approach 
of. fo zero as O " “ (z 1). 
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w we shall determine the multiplicative constant 
of the first eigenfunction of the longitudinal curvature 
problem. Again, the equation (5.62) is solved subject to the 
normalizing^ condition P^'(o) = 1. The first, eigenvalue is 
y, = 1.092. Equations (5.59), (5.60) and (5.61) 

are numerically integrated to yield the following solutions; 

^oo (0) = .4371, M^(0) = 56.416, p"(o) = -1.2344. (5.82) 

Using (5.57a) and (5.58) along with the values obtained above, 
tno left hand side of the equation (5.56a) is expressed in the 
following form for = 1/2; 

B* f"(o~, 0) .077cr~”’ + B^(2)~^/^ -1.546 + 0(o-“^). (5.85) 

The Eulerized Gortler series of the total skin friction due 
to longitudinal curvature *^2^^ as obtained in Chapter 4, is 
given by, 

= -1.352 + 2.014 Z -0.3750 Z^ - 0.1302 7? 

- 0.0600 Z"^. (5.84) 

We determine B^ by fitting (5.85) to (5.84). Since the proce- 
dure adopted is similar to the one outlined for the displacement 
speed problem, no details are given here. The value of the 
coefficient B^ is found to be -0.22. Substituting the value of 
B^ in (5.83), we obtain the following asymptotic series for 
the contribution of the longitudinal curvature to the skin 
friction for large values of o" ; 

B* F‘’^^^(cr7 0) .077 0"'"’’ - .039 + 0(cr~^), (5.£ 

The above series expansion shows that for large values of 
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curvature 

the contribution becomes aero iheee u ‘ 

With the earlier results bt • consistent 

shown in u ■ Eulerisea Gortler series 

snown in Figure 4,18. 


Combining 
asymptotic series 
thus, 


(5.73), (5.81) and ( 5 . 85 ), we obtain the 
for the total skin friction coefficient, 


Cf = 0.92768 + 0.11598 cr-"^ .. o .026 (r-‘^*546 

■n-1/2 / c + • . . 

tE (-0.017 0--5 ^ 0.077 cr -I ^ 0.289 <r-1-5 

- .039 o -"'**546 . 

( 5 . 86 ) 

The above asymptotic form of shows that the total skin 

u for - \f 2 approaches a value egual to 0.92768 as 
cr- ~ (a - 1 ). this behaviour was exhibited by the Improved 
Eulerlsed Gortler series obtadued to Chapter 4 and shown to 
figure 4.22b. Thus, the asymptotic series (5.86) provides, to 
some sense, a proof of convergence of the Bulerized GSrtler 
series solutions obtained In Chapter 4 and consequently establish 
the validity of the latter for the complete range of c . 


We, thus, conclude that the power series solution obtained 
for the nose region with the range of validity extended by employ, 
lug the Euler transformation, and the asymptotic solution for 
large streamwlse distance from the nose region provide means 
for assessing the accuracy of the corresponding results obtained 
via numerical or approximate methods having unknown and/or 
racer tain error bo rads. 
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Fig. 3.1 -Solutions to longitudinal curvature nnomentum transfer problem (3.14) 
En,n=1-4 denote the first four eigenvalues of the operator (3.16) 
for ,^ = 0.5 and Di , the first eigenvalue for ^=1.0 . 













recovery factor due to transverse curvature 






= 0.0,Pr=3.0 



Effect of displacement speed on heat transfer 
prescribed wall temperature. 





Fig. 3 . 9 - Second - order recovery factor due to displacement speed 




prescribed wall temperature . 




Fig. 3.11- Second “ order recovery factor due to vorticity interaction 



- order recovery factor du 
ng stream: Liquids. 
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Fig. 3.15- Second - order recovery factor due to displacement speed : Gases. 




Fig.4;i_A blunted wedge configuration with conformal mapping used 
to obtain the inviscid flow velocity. 
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Fig. 4.2c-First -order velocity profile on a blunted wedge 

I = 0 . 8 ) 





Fig. 4.4 -First -order skin friction distribution on a blunted wedge. 






Fig. 4.6 _ First - orders heat transfer on a blunted wedae 








Fig. 4.8c -Second - order velocity profile on a blunted wedge ( ^0^= 0.8 ) 
due to displacement effect . 





Displacement speed contribution to skin friction distribution 
on a blunted wedge . 
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Fig. 4.12c - Second - order temperature profile on a blunted wedge (^=0.8) 
due to displacement effect . °° 




Fig. 4.14 >. Displacement speed contribution to heat transfer distribution 
on a blunted wedge . 





Fig. 4.1 5b- Second -order velocity profile on a blunted wedge ( =0.5) 
due to longitudinal curvature °° 




Fig, 4.15 c -Second - order velocity profile on a blunted wedge (^08) 
due to longitudinal curvature. 









Fig. 4.17 -Contribution of longitudinal curvature to displacement thickness 
on blunted wedge. 










Fig. 4.19 c _ Second -order temperature profile on a blunted wedge (^=0.8) 
due to lonaitudinal curvature . °° 
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Fig. 4.20a- Contribution of longitudinal curvature to heat transfer 
on a blunted wedqe. 
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Fig. 4.22a- Skin friciion disiribution on a blunted wedge . 









Fig. 4.24 ^Comparison of present results of skin friction on a parabola 
(A=0) with exact solution .of the Navier- Stokes equations. 
















ions of equation ( 5.30) for 









